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If you fall into a real astronomical black hole (choosing a supermassive black hole, to make sure that the tidal forces don't 
get you first), then you will probably meet your fate not at a central singularity, but rather in the exponentially growing, 
■relativistic counter-streaming instability at the inner horizon first pointed out by Poisson & Israel (1990), who called it 
mass inflation. The chief purpose of this paper is to present a clear exposition of the physical cause and consequence 
of inflation in spherical, charged black holes. Inflation acts like a particle accelerator in that it accelerates cold ingoing 
and outgoing streams through each other to prodigiously high energies. Inflation feeds on itself: the acceleration is 
powered by the gravity produced by the streaming energy. The paper: (1) uses physical arguments to develop simple 
approximations that follow the evolution of inflation from ignition, through inflation itself, to collapse; (2) confirms 
that the simple approximations capture accurately the results of fully nonlinear one- and two-fluid self-similar models; 
1(3) demonstrates that, counter-intuitively, the smaller the accretion rate, the more rapidly inflation exponentiates; (4) 
shows that in single perfect-fluid models, inflation occurs only if the sound speed equals the speed of light, supporting 
the physical idea that inflation in single fluids is driven by relativistic counter-streaming of waves; (5) shows that what 
happens during inflation up to the Planck curvature depends not on the distant past or future, but rather on events 
happening only a few hundred black hole crossing times into the past or future; (6) shows that, if quantum gravity does 
not intervene, then the generic end result of inflation is not a general relativistic null singularity, but rather a spacelike 
singularity at zero radius. 
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1. Introduction 



It has been known for decades that something strange 
and dramatic must happen near the inner horizon of a re- 
alistic black hole. Penrose (1968) [l| page 222 first pointed 
out that a person passing through the Cauchy horizon (the 
outgoing inner horizon) of an empty charged (Reissner- 
Nordstrom) or rotating (Kerr-Newman) black hole will see 
the outside universe infinitely blueshifted. Empty here 
means empty of charge, matter, or radiation (but not of 
a static electric field) except at the singularity. Subse- 
quent perturbation theory investigations of the Reissner- 
Nordstrom geometry, starting with [2] and culminating 
with , showed that perturbations from the outside world 
would amplify to an infinite fiux of energy on the Cauchy 
horizon, a result that was interpreted as indicating that 
the Reissner-Nordstrom geometry was unstable. An early 
attempt to go beyond the linear regime is fi] . 

The full nonlinear nature of the instability near the 
inner horizon was eventually clarified in a seminal paper 
by Poisson & Israel (1990) ,5, 6]. Poisson & Israel argued 
that, if ingoing and outgoing streams are simultaneously 
present just above the inner horizon of a spherical charged 
black hole, then relativistic counter-streaming between the 
ingoing and outgoing streams will lead to an exponentially 
growing instability which they dubbed "mass inflation." 
During mass inflation, the interior, or Misner-Sharpe 0, 
mass, a gauge-invariant scalar quantity, exponentiates to 
huge values. Other gauge-invariant measures, such as the 
proper density, the proper radial pressure, and the Weyl 
scalar, exponentiate along with the interior mass. The 
phenomenon of mass inflation has been confirmed analyti- 
callyand numerically in many papers [1, 0, [lO, 11, 12, 13, 

Real black holes probably have very little charge, thanks 
to the huge charge-to-mass ratio of individual protons and 
electrons, e/m^ « 10^^, where e is the dimensionless charge 
of the proton or electron, the square root of the fine- 
structure constant, and nip is the proton mass in units 



of the Planck mass. A charged black hole would quickly 
attract opposite charge and almost neutralize, though the 
black hole would probably retain a small residual positive 
charge because protons are more massive than electrons, 
therefore more easily able to overcome a Coulomb barrier. 
Whereas real black holes have little charge, they probably 
do rotate (and rotate rapidly). The interior structure of 
a charged black hole resembles that of a rotating black 
hole in that the negative pressure (tension) of the elec- 
tric field produces a gravitational repulsion analogous to 
that produced by the centrifugal force in a rotating black 
hole. As a result, charged spherical black holes have inner 
horizons like rotating black holes. It is thus common in 
modeling black hole interiors to take charge as a surrogate 
for angular momentum, e.g. [^, H . 17 1, and the present 



paper follows this route. Further work on the interiors of 
rotating black holes is desirable, but the pro blem is hard, 
and investigations remain few 3, 1^ 2^, 21, 22] (see also 
the review j23i]). 

Much of the literature on mass inflation, starting with 
has focused on the situation where a black hole col- 
lapses, and then remains isolated. In this case, perturba- 
tions scattering off the curvature of the collapsing black 
hole are expected to lead to a decaying Price tail [24, 25l 
of outgoing gravitational radiation, as originally proposed 
by [ll- This outgoing radiation has usually been modeled 
with a massless scalar field_Eij23j28', '2?, 'so!^^ 

[HI, 0, S H m [H, S, lillialil iM, 39, io^iirEl lig, 

which supports waves moving at the speed of light, and 
is thus supposed to mimic gravitational radiation while 
retaining the simplifying advantage of spherical symmetry. 

Real astronomical black holes are however never iso- 
lated. Supermassive black holes at the centers of galax- 
ies are thought to build up their observed ^ 10^-10^ Mq 
masses 43|, |4J, |45| by gradual, albeit sporadic, accretion 
over the age of the Universe. Astronomical black holes, 
whether supermassive or stellar-sized, will continue to ac- 
crete some mixture of baryons and dark matter from their 
galactic environments. Even if nothing else is present, 
a black hole will accrete cosmic microwave background 
photons. 

Accretion provides an ongoing source of ingoing and 
outgoing matter inside a real black hole. Matter free- 
falling from outside the horizon is necessarily initially in- 
going, but it may become outgoing inside the horizon if 
it has enough charge (if the black hole is charged) or 
enough angular momentum (if the black hole is rotating) . 
It seems likely that outgoing radiation produced as a result 
of accretion will soon overwhelm the rapidly decaying Price 
tail of outgoing radiation produced by the initial collapse 
event. This is especially true in a supermassive black hole, 
whose mass acquired by accretion greatly exceeds the mass 
of the original black hole formed by stellar collapse. 

Occasional events of high accretion, most notably for 
example from the merger of two black holes, will re-energize 
a Price tail of outgoing radiation. In a general situation 
of variable accretion, outgoing radiation will be generated 
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by some combination of ring-down from peak accretion 
events, and steady accretion. 

The purpose of this paper is to clarify the physics of the 
mass inflation instability inside spherical, charged black 
holes. As originally proposed by [1], mass inflation is fueled 
by counter-streaming ingoing and outgoing streams. We 
take the point of view that ingoing and outgoing streams 
are being generated continuously by accretion. We treat 
the rates of generation of ingoing and outgoing streams 
as free parameters, and investigate how mass inflation 
depends on these rates. 

In this paper we assume that ingoing and outgoing 
streams are non-interacting. At typically low astronomical 
accretion rates, the assumption of non-interacting streams 
is likely to be a good one, at least in the realm of "known 
physics." However, the assumption of non- interacting streams 
is likely to break down at high accretion rates, or at center- 
of-mass collision energies exceeding the Planck mass. In- 
teraction between ingoing and outgoing streams is likely 
to have interesting consequences. We hope to explore the 
situation of interacting streams in a subsequent paper. 

The strategy of this paper is to start from qualitative 
and approximate arguments, and to proceed to illustrative 
models that support the qualitative conclusions. To put 
the paper into context, ^ presents the example of the 
supermassive black hole at the center of our Galaxy. To 
start the qualitative argument, identifies the two Ein- 
stein equations (|8]) that are at the heart of the mechanism 
of inflation. Next, 2] presents the simplest possible exam- 
ple, that of two symmetrically equal ingoing and outgoing 
streams, each neutral, pressureless, and freely-falling. The 
assumption of symmetrically equal ingoing and outgoing 
streams is equivalent, i i4.4l to the stationary approxima- 
tion d/dt = 0, previously introduced by 46|, |47|, |33|, who 
call it the homogeneous approximation because the time 
direction t is spacelike inside the horizon (we follow the 
convention of [48, p. 203] in referring to time translation 
symmetry as stationary, even when the time direction is 
spacelike rather than timelike). The assumption of equal 
streams is not realistic, so 311 goes on to generalize to the 
case of two unequal streams. 

The paper then moves on from approximations to self- 
consistent, self-similar models of accreting, charged, spher- 
ical black holes [H, Section [6] considers two-stream 
baryon-plus-dark-matter models, while Sj7]considers single- 
fluid models, including the case of a perfect fluid with an 
ultra-hard equation of state, which can be considered as a 
model of a massless scalar field 

All of this would be irrelevant if what happens during 
inflation depends on events that occur in the far future 
(as suggested by the Penrose diagram of the Reissner- 
Nordstrom geometry), so SJ5] discusses that issue. 

Section[2]checks whether the assumption of non- interacting 
streams is valid. 

To conclude, SJTO] summarizes the findings of the paper, 
and fJTT] appends some final remarks. 



2. An astronomical example 

It is useful to put the arguments of subsequent sections 
into context by considering a particular example: the ~ 
4 X 10^ MfD supermassive black hole at the center of our 
In this section, results from subsequent 



Galaxy [51|, |52| 
sections will simply be quoted without justification. 

Mass inflation takes place at (or rather just above) the 
inner horizon. Unless the black hole is rotating unusually 
slowly, the inner horizon of the Milky Way black hole will 
be at a radius r_ of 



GM/c^ - 10^ km , 



(1) 



or several times the radius of the Sun. 

The characteristic dimensionless accretion rate /i of the 
black hole, which can be thought of as the velocity with 
which the horizon is expanding in units of the speed of 
light, is expected to be of the order of the light crossing 
time of the black hole, about 1 minute, divided by the 
accretion timescale, about the age of the Universe, 



GM/c^ 



1 min 



age of Universe 10^" yr 



10" 



(2) 



which is tiny. 

Mass inflation is driven by relativistic counter-streaming 
between ingoing and outgoing streams. During inflation, 
the counter-streaming velocity u (meaning jv) exponenti- 
ates with an e-folding lengthscale of about /ir_ (eq. (fT9)) 
or IHI)), which is 



fir- 



10- 



(3) 



Thus as the counter-streaming ingoing and outgoing streams 
drop 1 millimeter in radius, the counter-streaming velocity 
increases by 1000 e-folds. At this point, still at a radius 
of ^ 10^ km, the Weyl curvature, center-of-mass energy 
density, and such-like quantities easily exceed the Planck 
scale. 

Presumably other physics of some kind intervenes, whether 
it be quantum gravity, or large collision cross-sections at 
super-Planck energies. But if not, if ordinary general rela- 
tivity continues to operate, and if the streams remain non- 
interacting, then the counter-streaming velocity would ex- 
ponentiate up to the absurdly huge value of (eq. (^2]) or (|T7| ) 



(4) 



before the inflationary period came to an end. 

During inflation, ingoing and outgoing streams both 
see each other highly blueshifted. An outgoing observer 
sees ingoing material accreted from the external Universe 
to the future of the time at which the observer jumped 
into the black hole, while an ingoing observer sees outgoing 
material accreted from the external Universe to the past of 
the time that they jumped in. Despite the huge blueshift, 
each stream sees only a modest amount of time go by on 
the other stream, approximately one black hole crossing 
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time per e-fold increase in blueshift (eq. ([52])). After 1000 
e-folds of blueshift, easily enough to pass the Planck scale, 
each stream has seen less than one day to have passed 
on the other stream. Thus what happens during inflation 
depends only on the immediate ftiture and past of the black 
hole, not on its distant future or past. 

From each of the ingoing and outgoing stream's own 
perspective, all this happens in a blink, approximately 
the light-crossing time across one e-folding length scale 
of 10"'^ mm, which is a few femtoseconds (10~^^s). Hit- 
ting the inner horizon is like hitting the pavement at the 
bottom of a skyscraper. You sailed down, but then it's all 
over in an instant. 

Although the center-of-mass energy density of the counter- 
streaming ingoing and outgoing streams is vast, the proper 
energy density of each stream individually remains modest, 
comparable to its energy density before inflation. Two 
cold streams being accelerated through each other up to 
highly relativistic velocities is characteristic of a particle 
accelerator. The situation may aptly be termed the "black 
hole particle accelerator" {^^■ 

3. The mechanism of mass inflation 

The purpose of this section is to elucidate the physics of 
mass inflation sufficiently well to provide a sound founda- 
tion for the approximations developed in the following two 
sections, fj4l ^-i^d ill The arguments and approximations 
are intended to be simple, but the reader should not be 
deceived into thinking that it is trivial or obvious to de- 
rive approximations that capture the behavior of inflation 
accurately. To date, the only simple a ppr oximation that 



M, M |33, 'TBI is 
which 1461 calls 



has been considered in the literature 
the stationary approximation d/dt = 0, 
the homogeneous approximation since the time direction 
is spacelike near the inner horizon. As shown ii4.4[ the 
stationary approximation is equivalent to the (unrealistic) 
case of equal ingoing and outgoing streams. 

This paper follows [i^, [3] in adopting the following 
line-element for a general spherically symmetric spacetime, 
in polar coordinates x'^ = {t, r, 9, (/)}: 



a 



-sin^eidf/)^) . (5) 



Here r is the circumferential radius, a gauge- invariant scalar, 
defined such that the proper circumference at radius r 
is 27rr. As detailed in Appendix [Appendix A[ the line- 
element ([5]) encodes not only a metric, but a complete 
orthonormal tetrad 7™ = {7t, 7r, 7e, 70}, that is, a lo- 
cally inertial frame, at each point of the spacetime. The 
coefficients /3t and Pr in the line-element ([5]) constitute 
the components of a tetrad-frame 4-vector, the radial 4- 
gradient 

/3m = dnir , (6) 

where = 7„i • d = e„/'(9/9x'' denotes the directed 
derivative along the 'jm tetrad axis (not to be confused 



with the coordinate derivative d/dx^^). Physically, the 
directed derivatives dt and dr are the proper time and 
radial derivatives measured by a person at rest in the 
tetrad frame. In particular, Pt = dtr equals the proper 
rate of change of the circumferential radius r measured by 
a person at rest in the tetrad frame (in another notation, 
this might be written dr/dr, where r is the proper time 
in the tetrad frame). The scalar length squared of /3m 
defines the interior, or Misner-Sharp [7], mass M{r), a 
gauge-invariant scalar, by 



2M 



1 = (3^ 



(7) 



For reference. Appendix [Appendix A| summarizes the var- 
ious entities, such as the tetrad connections and Einstein 
tensor, derived from the line-element (jS]). 

Spherically symmetric spacetimes are described by four 
independent Einstein equations. Two of the four serve 
to enforce conservation of energy-momentum, so if the 
equations governing the energy-momentum of the system 
are arranged to satisfy conservation of energy-momentum, 
as they should, then any two of the Einstein equations 
can be discarded as redundant. For the two remaining 
Einstein equations, it is most insightful to take the expres- 
sions (|A.10c|) and (|A.10b|) for the components C'' and G*'' 
of the Einstein tensor, which yield the following Einstein 
equations for the time evolution of jSt and Pr'- 



M 



A'Krp , 



(8a) 
(8b) 



where Dm denotes the tetrad-frame covariant derivative, 
/ = T*'' is the energy flux, and p = T^^ is the radial 
pressure measured in the tetrad frame. Equations (|S]) are 
valid in any (radially-moving) tetrad frame, and for any 
form of energy-momentum; the only requirement is that 
the spacetime be spherically symmetric. 

The two Einstein equations © are at the heart of mass 
inflation. Equation ([5a|) relates the radial acceleration on 
the left hand side to the gravitational force on the right 
hand side. The left hand side of equation ([5a|) . Dtfit, 
is the proper acceleration of the circumferential radius r 
measured by an observer who is in free-fall and instanta- 
neously at rest in the tetrad frame. The right hand side of 
equation (jSal) is the gravitational force, which consists of 
two terms, the (apparently) Newtonian gravitational force 
— M/r^, and an additional term —Anrp proportional to the 
radial pressure p. 

As will be seen below, it is the second, pressure term, 
in the gravitational force on the right hand side of equa- 
tion ((8aj) that is the source of the fun, first initiating mass 
inflation, and then driving it exponentially. 

Like the second half of a vaudeville act, the second 
Einstein equation (j8b[) also plays an indispensible role. 
The quantity /3r = drV on the left hand side is the proper 
radial gradient of the circumferential radius r measured 
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by a person at rest in the tetrad frame. The sign of /3r 
determines which way an observer at rest in the tetrad 
frame thinks is "outwards," the direction of larger circum- 
ferential radius r. A positive /3r means that the observer 
thinks the outward direction points away from the black 
hole, while a negative l3r means that the observer thinks 
the outward direction points towards from the black hole. 
Outside the outer horizon f3r is necessarily positive, be- 
cause must be spacelike there. But inside the horizon 

may be either positive or negative. A tetrad frame can 
be defined as "ingoing" if the proper radial gradient Pr is 
positive, and "outgoing" if is negative. In the Reissner- 
Nordstrom geometry, ingoing geodesies have positive en- 
ergy, and outgoing geodesies have negative energy. How- 
ever, the present definition of ingoing or outgoing based on 
the sign of is general - there is no need for a timelike 
Killing vector such as would be necessary to define the 
(conserved) energy of a geodesic. 

Equation (j8b[) shows that the proper rate of change 
DtPr in the radial gradient /3r measured by an observer 
who is in free-fall and instantaneously at rest in the tetrad 
frame is proportional to the radial energy flux / in that 
frame. But ingoing observers tend to see energy fiux point- 
ing away from the black hole, while outgoing observers 
tend to see energy flux pointing towards the black hole. 
Thus the change in j3r tends to be in the same direction 
as Pr, amplifying /3r whatever its sign. 

3.1. Reissner- Nordstrom phase 

Figure [T] illustrates how the two Einstein equations dSJ) 
produce the three phases of mass inflation inside a charged 
spherical black hole. 

During the initial phase, illustrated in the top panel of 
Figure [TJ the spacetime geometry is well-approximated by 
the vacuum (more correctly, electrovac) Reissner- Nordstrom 
geometry. During this phase the radial energy flux / 
is effectively zero, so /3r remains constant, according to 
equation (|8b|) . The change in the infall velocity /3t, equa- 
tion ([5a|) . depends on the competition between the Newto- 
nian gravitational force — M/r^, which is always attractive 
(tending to make the infall velocity j3t more negative), 
and the gravitational force —Airrp sourced by the radial 
pressure p. In the Reissner- Nordstrom geometry, the static 
electric field produces a negative radial pressure, or ten- 
sion, p = — (5^/(87rr^), which produces a gravitational 
repulsion —Airrp = Q'^/(2r^). At some point (depending 
on the charge-to-mass ratio Q/M) inside the outer hori- 
zon, the gravitational repulsion produced by the tension 
of the electric field exceeds the attraction produced by 
the interior mass M, so that the infall velocity Pt slows 
down. This regime, where the (negative) infall velocity /3t 
is slowing down (becoming less negative) , while remains 
constant, is illustrated in the top panel of Figure [T] 

If the initial Reissner-Nordstrom phase were to con- 
tinue, then the radial 4-gradient /3m would become light- 
like. In the Reissner-Nordstrom geometry this does in 




Fi gure 1: Spacetime diagrams of the tetrad- frame 4- vector /?m, 
equation ||6ll, illustrating qualitatively the three successive phases 
of mass inflation: 1. (top) the Reissner-Nordstrom phase, where 
inflation ignites; 2. (middle) the inflationary phase itself; and 
3. (bottom) the collapse phase, where inflation comes to an end. 
In each diagram, the arrowed lines illustrate two representative 
examples of the 4-vector {/3t,/3r}, one ingoing and one outgoing, 
while the double- arrowed lines illustrate the rate of change of these 
4- vectors implied by Einstein's equations ||8ll. Inside the horizon of 
a black hole, all locally inertial frames necessarily fall inward, so the 
radial velocity /3t = dtr is always negative. A locally inertial frame is 
ingoing or outgoing depending on whether the proper radial gradient 
/3r = drr measured in that frame is positive or negative. 



fact happen, and where it happens defines the inner hori- 
zon. The problem with this is that the lightlike 4-vector 
f3m points in one direction for ingoing frames, and in the 
opposite direction for outgoing frames. If (3m becomes 
lightlike, then ingoing and outgoing frames are streaming 
through each other at the speed of light. This is the infinite 
blucshift at the inner horizon first pointed out by [Ij. 

If there were no matter present, or if there were only 
one stream of matter, either ingoing or outgoing but not 
both, then /3„i could indeed become lightlike. But if both 
ingoing and outgoing matter are present, even in the tiniest 
amount, then it is physically impossible for the ingoing and 
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outgoing frames to stream through each other at the speed 
of Ught. 

If both ingoing and outgoing streams are present, then 
as they race through each other ever faster, they generate 
a radial pressure p, and an energy flux /, which begin to 
take over as the main source on the right hand side of 
the Einstein equations ([8|). This is how mass inflation is 
ignited. 

3.2. Inflationary phase 

The infalling matter now enters the second, mass infla- 
tionary phase, illustrated in the middle panel of Figure [1] 

During this phase, the gravitational force on the right 
hand side of the Einstein equation ([5a|) is dominated by the 
pressure p produced by the counter-streaming ingoing and 
outgoing matter. The mass M is completely sub-dominant 
during this phase (in this respect, the designation "mass 
inflation" is misleading, since although the mass inflates, it 
does not drive inflation). The counter-streaming pressure 
p is positive, and so accelerates the infall velocity (3t (makes 
it more negative). At the same time, the radial gradient 

is being driven by the energy flux /, equation (I8bp . 
For typically low accretion rates, the streams are cold, 
in the sense that the streaming energy density greatly 
exceeds the thermal energy density, even if the accreted 
material is at relativistic temperatures. This follows from 
the fact that for mass inflation to begin, the gravitational 
force produced by the counter-streaming pressure p must 
become comparable to that produced by the mass M, 
which for streams of low proper density requires a hyper- 
relativistic streaming velocity. For a cold stream of proper 
density p moving at 4-velocity u™ = {u*,u^, 0,0}, the 
streaming energy flux would be / ^ pu^vT, while the 
streaming pressure would he p ^ p(u'")^. Thus their ratio 
f Ip ^ /u^ is slightly greater than one. It follows that, 
as illustrated in the middle panel of Figure [TJ the change 
in Pr slightly exceeds the change in Pt , which drives the 4- 
vector f3m, already nearly lightlike, to be even more nearly 
lightlike. This is mass inflation. 

Inflation feeds on itself. The radial pressure p and 
energy flux / generated by the counter-streaming ingo- 
ing and outgoing streams increase the gravitational force. 
But, as illustrated in the middle panel of Figure \T\ the 
gravitational force acts in opposite directions for ingoing 
and outgoing streams, tending to accelerate the streams 
faster through each other. As ^16.] put it, the gravitational 
force is always inwards, meaning in the direction of smaller 
radius, but the inward direction is towards the black hole 
for ingoing streams, and away from the black hole for 
outgoing streams. 

The feedback loop in which the streaming pressure and 
flux increase the gravitational force, which accelerates the 
streams faster through each other, which increases the 
streaming pressure and flux, is what drives mass inflation. 
Inflation produces an exponential growth in the streaming 
energy, and along with it the interior mass, and the Weyl 
curvature. 



3.3. Collapse phase 

It might seem that inflation is locked into an exponen- 
tial growth from which there is no exit. But the Einstein 
equations ([8]) have one more trick up their sleave. 

For the counter-streaming velocity to continue to in- 
crease requires that the change in /3r from equation (jSbp 
continues to exceed the change in /3t from equation (j8ap . 
This remains true as long as the counter-streaming pres- 
sure p and energy flux / continue to dominate the source 
on the right hand side of the equations. But the mass 
term —M/r^ also makes a contribution to the change in fit. 



equation ((Sa)) . As will be seen in the examples of the next 
two sections, §Sj31a'nd[51 at least in the case of pressureless 
streams the mass term exponentiates slightly faster than 
the pressure term. At a certain point, the additional ac- 
celeration produced by the mass means that the combined 
gravitational force M/r^ + Airrp exceeds inrf. Once this 
happens, the 4- vector /Sm, instead of being driven to be- 
coming more lightlike, starts to become less lightlike. That 
is, the counter-streaming velocity starts to slow. At that 
point inflation ceases, and the streams quickly collapse to 
zero radius. 

It is ironic that it is the increase of mass that brings 
mass inflation to an end. Not only does mass not drive 
mass inflation, but as soon as mass begins to contribute 
significantly to the gravitational force, it brings mass in- 
flation to an end. 

4. Two equal streams 

The purpose of this section is to present the simplest 
possible quantitative example of how the counter-streaming 
instability described qualitatively in the previous section 
works. The model is that of two equal ingoing and out- 
going streams of neutral, pressureless matter. As shown 
in ji4.4[ the approximation is equivalent to the stationary 
(homogeneous) approximation first considered by 46|, |47 , 
133| . The assumed symmetry between ingoing and outgoing 
streams is not realistic, but has the virtue of simplicity. A 
more realistic but still idealized model with two unequal 
streams is solved in the next section, ^ Self-consistent 
models are presented in later sections. 

For pedagogy, the details of the solution for two equal 
streams are left as a problem for the reader, Appendix |Appendix B 
This problem appeared on the final exam of the first au- 
thor's graduate course on general relativity in Spring 2008. 
The problem includes sufficient hints that the reader should 
be able to solve it. In any case, elements of the solution 
are given below. 

As seen in the previous section, ^ the key object 
is the radial 4- gradient (3m, equation The radial 4- 
gradient contains two degrees of freedom, its time and 
radial components I3t and f3r- It is convenient to recast 
these two degrees of freedom in terms of a magnitude 
P = I /3m I and the radial component u = li'' of a radial 
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4-velocity u" 



Table 1: Phases of evolution. 



l3,n = {/3uPr,0,0} = ~l3Wl + u\u,0,0} . (9) 

The magnitude f3, which is positive, is related to the in- 
terior mass M by /S^ = 2M/r - 1, equation The 4- 
velocity u™ = + u^, u, 0, 0} is the 4-velocity of the 
stream relative to the no-going frame at the boundary 
between ingoing and outgoing, where Pr = 0. The sign 
oi u = is negative for ingoing, positive for outgoing. 

The central equation that emerges from the analysis in 
Appendix [Appendix B| is [the following equation repeats 
eq. HJl)] 



dln/3 



(10) 



dhiu X — P"^ + fi + /iu^ 

The dimensionless parameters A and fi in equation (|10p 
are well-approximated as constant during mass inflation. 
The dimensionless parameter A is 



A = gVr2 - 1 



(11) 



where Q is the interior charge of the black hole and r is the 
radius. It is a feature of inflation that during it the radius 
r hardly changes, a fact confirmed below, equation (jl9p . 
The radius may be approximated by the radius of the inner 
horizon r « r_ — M ~ ^ A'P — of the parent Reissner- 
Nordstrom black hole. The parameter A is typically a 
number of order unity; for example, if Q/M = 0.8, then 
A = 3. The dimensionless parameter in equation (jlOp is 



(12) 



proportional to the product of the radius r squared and 
the proper mass-energy density p of either of the ingo- 
ing or outgoing streams. Perhaps surprisingly, the proper 
density p of either stream in its own frame remains almost 
constant during inflation. What exponentiates during in- 
flation is the counter-streaming velocity, and along with 
it the center-of-mass energy density; by comparison the 
energy density of each stream changes little, because its 
volume changes little. The fact that a volume element 
remains little distorted during inflation even though the 
tidal force, as measured by the Weyl curvature scalar, 
exponentiates to huge values was first pointed out by Q. 
The physical reason for the small tidal distortion despite 
the huge tidal force is that the proper time over which the 
force operates is tiny. The parameter p is proportional to 
the accretion rate of the black hole, and for simplicity p 
is referred to loosely in this paper as the "accretion rate." 
As seen in ^J2j astronomical accretion rates are typically 
small (tiny), so p is small (tiny). 

Equation (jlOp is a one-dimensional ordinary differential 
equation, whose solution defines curves in the plane. 
The behavior depends on the relative magnitudes of three 
terms, A, pu^ , and (a fourth term p is never important) 
that appear in the numerator and denominator on the 
right hand side of equation (flUl) . The three phases of 
evolution described in fJ3] occur as each of the three terms 



Phase 


Gravity 


Dominant term 


Reissner- Nordstrom 


charge 


A 


Inflation 


streaming 


pv?' 


Collapse 


mass 





dominates, Table [T] As listed in the middle column of 
Table [TJ physically, the dominance of each term expresses 
the dominance of three sources of the gravitational force 
on the right hand side of the Einstein equations ([5]): A 
dominates if the negative pressure of the electric field dom- 
inates; pv? dominates if the counter-streaming pressure 
and flux dominate; and dominates if the interior mass 
M dominates. The three phases are considered in the next 
three subsections. 

^.1. Reissner- Nordstrom phase 

Because the accretion rate p is small, initially the stream- 
ing pv?' term (and the p term) in equation pop can be 
ignored. During this initial, Reissner-Nordstrom, phase, 
equation ([TU)) simplifies to 



d\nl3 _ + 
dluM ^ A - /32 



= -1 



(13) 



The initial values of /3 and u are set by the boundary condi- 
tions of accretion, and are generically of order unity; their 
precise initial values are not important to the argument. 
The solution of equation is 



u 



(14) 



where C is a constant of order unity, set by the boundary 
conditions. This Reissner-Nordstrom phase of evolution is 
labeled "RN" in the fi-u diagram of Figure [2J 

In the Reissner-Nordstrom geometry, (3 would decrease 
to zero at the inner horizon, while the velocity u would 
diverge, in accordance with equation (ITi)) . This reflects 
the tendency of the 4-vector to become lightlike at the 
inner horizon, as described in SJ2] and illustrated in the top 
panel of Figure [TJ 

As the velocity u becomes larger on the approach to 
the inner horizon, the term pv?' in equation (1101) , which is 
proportional to the streaming pressure in the no-going (= 
center-of-mass, here) frame, eventually starts to become 
appreciable, however tiny the accretion rate p may be 
initially. The Reissner-Nordstrom phase ends when pvP' « 
A, which corresponds to the streaming pressure and flux 
taking over as the dominant source of gravity on the right 
hand side of the Einstein equations ([5]). 

4. 2. Inflationary phase 

During the Reissner-Nordstrom phase, /? is driven to 
a small value, while the velocity u increases to the point 
where the streaming pressure begins to take over. Since 
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Figure 2: The variation of /3 and u predicted by the simple model 
of two equal, pressureless, ingoing and outgoing streams. The 
parameters are A = 3 and fi = 0.1, and /3 and u are both initially 1. 



the terms in equation (|10p are now small, the equation 
may be approximated by 

- A + fiu"^ 



dln/3 
dlnu 



A + /iu^ 



(15) 



The solution of equation dTSl) that initially matches the 
Reissner-Nordstrom solution (ITil) is 



Xu 



(16) 



Equation (|15p shows that f3 reaches a minimum at = 
A/^, where, according to the solution (|T6| . 

/3-2|C|v/^. (17) 

The solution then enters the inflationary phase, where /3 
increases in proportion to u 



for li^ >• X/fi 



(18) 



This phase is labeled "inflation" in the /3-u diagram of 
Figure [2] 

That something dramatic is happening can be seen 
from the equations (|B.7|) relating the change in /3 and u to 
the change in the radius r. During the inflationary phase 
where the ^u^ term dominates, the inflationary growth 
rate from equations (jB.7|) is 



din (3 d In It 
d In r d In r 



A2 



2C2/1 



(19) 



Since fi is small (tiny), these rates are large (huge). That 
is, /3 and u increase dramatically while the radius r de- 
creases only a small amount. 

Equation (IT^ shows that the smaller the accretion rate 
II, the larger the inflationary growth rate. This counter- 
intuitive feature is confirmed in the self-consistent models 
described in §^JS]and[71 



4-3. Collapse phase 

During the inflationary phase the streaming term fj,u^ 
in equation (fTUl) is driven to a huge value, dwarfing the 
A term. But /3 has also been increasing during inflation, 
and the mass term in equation ^TU\i begins to make its 
presence felt. With the fxu^ and terms included but the 
A terms dropped, equation (1101) is approximately 



dln/3 _ + 
dlnu — + /iu^ 



(20) 



The right hand side of equation ([20]) is slightly greater 
than 1, which shows that /? grows slightly faster than u. 
Eventually, the term grows to the point that it equals 
IJ,u^. At this point the denominator of the right hand side 
of equation ([20)1 is zero. The velocity u then begins to 
decrease, while /3 continues to increase. This signals the 
end of mass inflation. 

Equation (j20[) has an exact solution given by the im- 
plicit equation 



R fx' 



exp 



in which the constant integration factor X'/ {Cii)e^ I '^'^'^ 
comes from matching to the inflationary solution (1181) at 
ff' = A, where equation (pHl) takes over from equation (fTS)) . 
Equation (|20l) shows that the velocity u reaches a maxi- 
mum when = ^u^, at which point equation (PT|) shows 
that 

A ( X' \\ 

Since the accretion rate /i is small, the maximum value 
of the velocity u is exponentially vast, the exponential 
of a large number, u ^ e^l^ . In realistic situations, the 
curvature and centcr-of-mass energy are likely to exceed 
the Planck scale before this vast value is attained, and 
quantum gravity will presumably intervene. 

However, if the accretion rate is large enough, as 
may happen during the first blush of gravitational col- 
lapse of a black hole, or during a black hole merger, then 
the maximum value (j22p may be attained before quantum 
gravity becomes important. It was argued in ^14. 21 that 
the radius r decreases only slightly during mass inflation. 
However, the number of e-folds that elapse in attaining 
the peak velocity ([^ is large enough that the radius 
drops appreciably. Specifically, in reaching the maximum 
velocity ([22)) the radius r decreases by approximately 



Alnr 



2CV 



A Inu : 



1 

2 ' 



(23) 



or a factor of e^^/^ « 0.6. 

After the velocity u peaks out, the matter collapses 
precipitously to zero radius. After the peak, equation (PT|) 
simplifies to 



y2 / y2 



(24) 



which is just a constant. This phase is labeled "collapse" 
in the (3-u diagram of Figure [21 During the collapse phase, 
the mass term dominates, and it follows from equa- 
tion (IB. Tap that /3 varies with radius r as 



din (3 
dlnr 



(25) 



This shows that the radius r is no longer stuck at its infla- 
tionary value. The proper radial infall velocity /3t « —f3u 
is exponentially vast, and the streams of matter plunge to 
zero radius in an instant. Since the interior mass M is 
related to /3 by 2M/r — 1 — /?^, and /3 is exponentially 
vast, equation (pS)) shows that 



dhiM 
dlnr 



^ 



(26) 



In other words, the interior mass M stops increasing. The 
value of M attained as the radius drops to zero is the 
exponentially huge value 



(27) 



where r_ , the radius of the inner horizon of the original 
parent Reissner-Nordstrom black hole, is the approximate 
value of the radius during inflation. 

4.4- Stationary approximation 

It might seem that in the limit of small accretion rates 
it would be reasonable to assume the stationary approxi- 
mation 

r) 

(28) 



since as seen from outside the black hole hardly changes 
as time goes by. However, as will now be shown, the 
stationary approximation (|28p can be applied consistently 
to mass inflation only in the (unrealistic) case of sym- 
metrically equal ingoing and outgoing streams, which is 
the case considered in this section. The approximation 
is termed the homogeneous approximation by 4^ 47. 33 1 
because the Killing vector associated with time translation 
symmetry is spacelike inside the black hole. We follow (isl , 
p. 203] in referring to the approximation as stationary, 
since it is just a single global symmetry of the spacetime. 
The conclusions of this section are consistent with those 
of [i^, 47, m, who emphasize that the stationary 



(homogeneous) approximation has the virtue of simplicity, 
and captures the general behavior of mass inflation, but 
does not predict accurately the behavior of more realistic 
models, even in the limit of vanishingly small accretion 
rates. 

The implications of the stationary approximation (1281) 
can be drawn from the following general identity, valid in 
arbitrary spherically symmetric spacetimes: 



A/3, {Gtt + Grr) - {lil + Pi) Gtr = 



2af3r dM 
'~ df 



(29) 



In the stationary approximation (I28p . the right hand side 
of equation ([^ vanishes, and equation ([^ then implies 
that in the center-of-mass frame, where Gtr — 0, one of 
the following three possibilities applies: 



either /3t = , 
or /3r = , 
or p + p = , 



(30) 



where p and p denote the proper center-of-mass density 
and radial pressure. Two of the three possibilities (|30)) . 
namely /3t = or p + p = 0, are incompatible with mass 
inflation. In the superluminally infalling inflationary zone, 
Prn is timelike, /3f > /3^, and (3t is strictly negative; indeed, 
— = 2M/r — 1 exponentiates to huge values during 
inflation, so (3t is hugely negative. Similarly, the proper 
center-of-mass density and radial pressure exponentiate to 
huge values during inflation, so p + p is hugely positive, 
not zero. 

This leaves the third possibility, that /3r — in the 
center-of-mass frame (note that af3r remains finite when 
Pr ~ 0)- The possibility corresponds to the case of sym- 
metrically equal ingoing and outgoing streams, where by 
symmetry the center-of-mass frame is the no-going frame, 
Pr = 0, at the boundary between ingoing and outgoing. As 
[Tij l put it, the stationary (homogeneous) approximation 
applied to the case of a massless scalar field implies that 
"the scalar field can be represented as a sum of two equal 
fluxes moving in opposite directions along the (spacelike) 
t axis." [15] also apply the homogeneous approximation 
to the case of single perfect fluids with either pressureless 
{p = 0) or relativistic {p/ p = 1/3) equations of state. How- 
ever, these single-fluid cases are irredeemably unrealistic, 
since they suppress the inflation that would occur if even 
the tiniest admixtures of ingoing and outgoing parts were 
admitted. 

The stationary approximation (j28p correctly reproduces 
the qualitative behavior of mass inflation, but does not 
capture accurately what happens in the more realistic case 
of unequal streams, fJS] If the stationary approximation (j28p 
were true, then geodesies would be characterized by a 
conserved energy (per unit mass) E', given by the covariant 
time component Uf of the coordinate-frame 4-velocity 
(the coordinate- frame 4-velocity v^^ is written with a Greek 
upsilon to avoid any possible confusion with the tetrad- 
frame 4-velocity u™, written with a Latin u). The ratio 
Eb /Ea of conserved energies of freely- falling streams h and 
a, expressed in the tetrad-frame of freely-falling stream a, 
is 



E^ 

Ea 



Vb£ 
Va,t 



Pa, 



Ua 



(31) 



In the situation considered in fJS] of two unequal streams 
accreting at constant rates, the ratio Eb/Ea of energies 
should be constant, and then equation (j3ip would imply 
that Ub/ua is constant. But in fact the inflationary solu- 
tion shows that Ub/ua is not constant. Rather it varies 



from U}j/ua = Ch/Ca in the Reissner-Nordstrom phase, 
equation (155)) . to u-^jua = Caiia/ (Cbiib) in the inflationary 
phase, equation (|1T|) . to 



Ub_ 

Ua 



CbfJ-b 



exp 



(32) 



in the collapse phase, equation (|44|) in the limit of small 
accretion rates and late times, Ua ^ 13. 

5. Two unequal streams 

The situation of symmetrically equal ingoing and out- 
going streams considered in the previous section, SjJl is 
straightforward to solve. However, the situation is not 
realistic, since in reality the mix of ingoing and outgoing 
streams depends on boundary conditions outside the outer 
horizon that are unlikely to produce equal streams. 

This section considers the more realistic situation of 
two unequal streams. The case remains analytically solv- 
able provided that, as in gl the streams are idealized as 
neutral and pressureless, and certain slowly- varying pa- 
rameters (A and /i) are again approximated as constants 
for each stream. The analytic solution for unequal streams 
can be compared to results from self-consistent models of 
accreting black holes, and it will be seen that the an- 
alytic solution captures accurately what happens in the 
self-consistent models. 

For two unequal streams, equation (jlOp generalizes to 
the two equations 



din/ 
dlnv 



A-/32 + 



in which a represents either of the ingoing or outgoing 
streams, and b represents the other stream. Equation (1331) 
is the ratio of equations (jC.7aP and (jC.7a[) derived in Ap- 
pendix Appendix C[ The accretion rate /it, of stream b in 
equation ([55)) is given by equation ([T^ . but subscripted. 



fib = 167rr pb . 



(34) 



Being a gauge- invariant scalar, (3 in equation (|55)) is the 
same for both streams. 

As in the case of equal streams considered in 21 the 
infalling matter goes through three phases of evolution. 
Table [TJ which are considered in turn in the next three 
subsections. 

5.1. Reissner-Nordstrom phase 

Realistically, the accretion rates pa of both ingoing and 
outgoing streams are likely to be small (even if one is much 
larger than the other), so the /ifjit^ terms on the right hand 
side of equation ()33l) can be ignored. This is the initial 
Reissner-Nordstrom phase, and during it equation p3p for 
each stream is well- approximated by 



dln/3 -A + /32 



dlnur, 



= -1 



(35) 



Equation (1351) solves to 



Ca 



(36) 



for each stream, where Ca is a constant, which will gen- 
erally be different for each stream. Note that Ca has 
the same sign as Ua, which is negative ingoing, positive 
outgoing. 

5.2. Inflationary phase 

As in the case of equal streams, /3 is driven to a small 
value during the Reissner-Nordstrom phase, while the ve- 
locities Ua of both streams increase, equation (155)) . Even- 
tually the fJ.b'U'l streaming terms on the right hand side 
of equation ([33]) become important, even though the ac- 
cretion rate /ib may be tiny. Since the mass terms 
in equation ([55)) are by that time small, equation 
becomes 

dln/3 —X + pbul 



dint 



A + iibui ^ 

where as before a represents either of the two streams, 
while b represents the other. Equation (1371) is a pair of cou- 
pled ordinary differential equations. The solution match- 
ing the Reissner-Nordstrom initial conditions ()36p is 

_ |Ca|(A + ^Jiaul) + [Ca{X - PaUD] ' + Xfib{2CbUa)^ 



2X\Ua\ 

(38) 

for either stream. Equation ([57)) shows that (3 goes through 
a minimum when /i^u^ = A, at which point /3 is, according 



(33) to the solution (1551) . 



f3 = 



\Ca\^/Jj^+ \Cb\^/JIb 

7x 



(39) 



which may be compared to the equal stream result ()17p . 



Given the solution ([58|. integration of equations (IC.7P 
yields the radius of each stream: 



(40) 



where the constant of proportionality may be set, for ex- 
ample, by the value of ra where /? goes through its min- 
imum. Note that the radii ra are (slightly) different for 
each stream, since each follows a different trajectory into 
the black hole. If Ua and Ub are expressed in terms of /3 
through equations ([38)) . then equation ([40)) gives the radius 
Ta of stream a as a function of 13 along its worldline. 

The solution then enters the inflationary phase, where 
13 and the velocities Ua and Ub increase in proportion to 
each other 

/3 = ^^t^ . (41) 
A 

The inflationary growth rate of each stream is, from equa- 
tions ([CJ)). 
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Figure 3: Absolute values |n| of velocities of outgoing baryonic (light 
blue) and ingoing dark matter (dark brown) streams as a function of 
/3 in a self-similar solution with charged, almost pressureless baryons 
and neutral, pressureless dark matter. The superposed dashed 
lines are the analytic approximation given by equation II38I I. The 
parameters of the model and of the fit are listed in Table [2] 



where as above a represents either of the two streams, and 
h the other. 

Figure [3] compares the analytic solution to a self- 
consistent self-similar solution of the kind described by 
[l6| and explored further in SjH For the model shown, the 
boundary conditions at the outer sonic point are: black 
hole charge-to-mass Q./M, — 0.8, mass accretion rate 
M, = 0.01, baryonic equation of state w = Pb/pb = 10^^, 
and proper dark-matter-to-baryon density Pd/pb — 0.001, 
as listed in Table[2j The model does not satisfy all the con- 
ditions of the analytic solution: the baryons are charged 
not neutral, and they are not quite pressureless. The 
baryons are charged for two reasons. First, in the self- 
similar solutions the charge of the black hole is created self- 
consistently by accretion, so at least one stream must be 
charged. Second, mass inflation requires the simultaneous 
presence of ingoing and outgoing streams. If both streams 
are accreted from outside the black hole, then one of the 
streams must be charged in order that, repelled by the 
charge of the black hole, the stream becomes outgoing 
inside the outer horizon. The baryons are not quite pres- 
sureless, having equation of state w = 10^^. The reason 
for this is that the boundary conditions of the model are 
set at a sonic point outside the black hole where the in- 
falling baryons accelerate from subsonic to supersonic. For 
pressureless baryons, it; = 0, the sonic point is at infinity. 
Giving w a small but finite value brings the sonic point to 
a large but finite radius. 

Figure [3] shows that, despite the fact that the model 
does not fulfill all the assumptions of the analytic solution, 
the fit is nevertheless very good (within the width of the 



line) just before and during inflation. Well before inflation, 
the baryons depart from the fit because they are charged, 
and therefore do not follow a geodesic. The parameters 
of the fit are given in Table [51 The parameters Ca of the 
analytic solution p8l) , namely Cb for the outgoing baryons 
and Cd for the ingoing dark matter, are fits to the Reissner- 
Nordstrom phase, equation (j36p . entering inflation. The 
values of the remaining parameters A, /if,, and pd of the 
analytic solution ([55)1 are given by equations (fTTI) and 
measured from the model at the onset of inflation, 
where /3 goes through its minimum. The parameters A, 
Pb^ and pd are perhaps not what might have been naively 
expected given the boundary conditions at the sonic point; 
a discussion of why they do in fact make sense is relegated 
to Appendix I Appendix D[ 

5.3. Collapse phase 

During inflation, the streaming terms PbU^ on the right 
hand side of equation p3p become huge, dwarfing the A 
terms. But infiation also drives /? to large values, eventu- 
ally to the point where the mass terms on the right hand 
side of equation (1551) can no longer be neglected compared 
to the pbul streaming terms. With the pbul and terms 
included, equation (j33l) becomes 



dln/3 /32 + 



Pbul 



din Ua 



-/32 + fl^ui 



(43) 



Equations (|43| constitute a pair of coupled ordinary dif- 
ferential equations, and they have an exact solution given 
by the implicit equations 



f3Ua 



A2 



Capa 



qt + Paul/p^ 



where qa are constants 



qb + XViClPa) 



+ Clpa 



9a /2 



CaPa 



ClPb 



(44) 



(45) 



and as before a represents either stream, and b the other. 
Note that qa + qb = 1, and that one of qa or g;, is < 
and the other is > 1. For small accretion rates /ia and 
Pb, the constants qa and qb are large and of opposite sign, 
so that the power law on the right hand side of equa- 
tion (I44p approximates an exponential. The overall factor 
in equation (j44p comes from matching to the inflationary 
solution (jlTl) at /S^ = A, where equation (|^51) takes over 
from equation (j37p . Integration of equations (|C.7p yields 
the radius r„ of each stream: 



\Ua\/P 



(46) 



where the constant of proportionality may be obtained 
by matching to equation (j40p with the inflationary solu- 
tion (HH) at = A. 
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Table 2: Model and fit parameters. 



Model parameters Fit parameters 
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Equation (^5)) shows that the velocity Ua reaches a 
= Pbul, at which point equation (I44p 



maximum when 
shows that 



.1/4 



^1/2 3/4 
L-a pa 



Qb 



9a /4 



(47) 



where Xa = —qb/{l + Qa) is the value of ^laU^/ fi"^ at the 
peak of Ua- 

After the velocity peaks, the matter collapses rapidly to 
zero radius. The interior mass, equation ([7]), is M = \0^ra 
to a good approximation, which according to equation ([^^ 
is proportional to fiua^ which from equation (|44)) tends to 
a constant at zero radius. In the typical case of small pa 
and pb^ the value of the interior mass M attained at zero 
radius is 



M ; 



■ exp 



\Hn{ClpalClpb) 



nClPa 



ClPb) 



(48) 



which generalizes the equal stream result ([27)) . 

Figure 0] compares the analytic solutions ((38)) and (|44)) . 
valid respectively during the RN+inflation and inflation+ 
coUapse phases, to a self-consistent self-similar model whose 
parameters are listed in Table [H Reaching the collapse 
phase in a model requires a large accretion rate /i, oth- 
erwise the peak velocity u ~ e^^^ and related quantities 
overflow numerically. A large accretion rate is more easily 
accomplished in models where the baryons have pressure, 
since then the sonic point is closer to the horizon. To 
achieve the desired large accretion rate, the model shown 
in Figure |4] has an almost relativistic baryonic equation 
of state w = 0.32 (a value of w slightly less than 1/3 
allows for an expected modest increase in the number of 
relativistic particles species as temperature increases) , and 
a large dark-matter-to-baryon density Pdj Pb = 0.1 at the 
sonic point. 

Figure H] shows that the approximate solutions provide 
a good fit to the behavior of the neutral, pressureless dark 
matter, and a satisfactory fit to the charged, relativistic 
baryons from just before inflation up to the end of infla- 
tion. As previously remarked for the model of Figure [31 
the fit to the baryons fails well before inflation because the 
baryons are charged, and therefore do not follow geodesies. 
Figure m also shows that the fit to the baryons fails as soon 
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Figure 4: Similar to Figure [S] illustrating all three phases of 
evolution in a self-similar solution similar to that in Figure |3] but 
in which the baryons, instead of being almost pressureless, have a 
relativistic equation of state w = Pb/Pb = 0.32. The superposed 
dashed lines are the analytic approximation II38II for the Reissner- 
Nordstr6m-|-inflation phases, and the analytic approximation II44II 
for the inflationary+coUapse phases. 



as the collapse phase begins, a fact that can be attributed 
to the pressure of the baryons, since at this point the 
Lorentz force on the baryons is completely subdominant. 
The parameters of the fit are listed in Table [H 

The failure in Figure U of the fit to baryons in the col- 
lapse phase is not terribly important. What is important is 
that the analytic approximation (j44p captures reliably the 
behavior during infiation, and it predicts approximately 
correctly when inflation comes to an end, equation (j47p . 
Once collapse begins, the ingoing and outgoing streams 
collapse to zero radius in a tiny instant. Collapse to zero 
radius is subject to t he non-radial Belinski-Khalatnikov- 
Lifshitz instabilities 



53l l54l . l55l |56|. Such instabilities. 



which are artificially suppressed as long as spherical sym- 
mmetry is assumed, and are not explored in this paper. 
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Figure 5: Similar to Figure [S] illustrating a self-similar solution 
in which one stream, the ingoing dark matter stream (upper dark 
brown line), is much smaller than the other, the outgoing baryonic 
stream (lower light blue line). The superposed dashed lines are the 
analytic approximation II38I I. 



5.4- Very unequal streams 

What happens if one of the two streams is very much 
smaller than the other? 

Figure [5] illustrates a self-consistent self-similar model 
in which the dark-matter-to-baryon density has been set 
to a low value Pd/pb = 10~^ at the sonic point. Figure [S] 
shows that the smaller dark matter stream undergoes pre- 
inflationary acceleration that brings its velocity Ud up to 
the point that CdP-dUd ~ CbfJ.bUb^ in accordance with equa- 
tion (|4ip. whereupon inflation ensues. The parameters of 
the model, and of the fit to it, are listed in Tabled 

According to equation (j48p , the maximum interior mass 
attained is of the order of 

M - e^/^''^ (49) 

where fj,a is the larger of the two accretion rates. 

6. Baryon + dark matter models 

The self-consistent models referred to in previous sec- 
tions are the general relativistic, self-similar, accreting, 
spherical, charged black hole solutions described by [49|, 
[16[, to which the reader is referred for more detail. In 
these models the black hole accretes two fluids: charged 
"baryons" with equation of state pt/pt = w] and neutral 
pressureless "dark matter." The baryons, being repelled 
by the charge of the black hole produced self-consistently 
from the accreted baryonic charge, become outgoing inside 
the horizon. The dark matter, being neutral, remains in- 
going. Relativistic counter-streaming between the baryons 
and the dark matter produces mass inflation. 




.05 .1 .2 .5 12 
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Figure 6: Interior mass M, equation ([TJ, as a function of baryonic 
radius rj, for models with dark-matter-to-baryon ratio Pdl Pb = 0.1 
at the sonic point, and three different accretion rates, M, = 0.003, 
0.01, 0.03. The smaller the accretion rate, the faster the interior 
mass M inflates. The superposed dashed lines are the analytic 
approximation II40I I for the Reissner-Nordstrom-f inflation phases, 
and the analytic approximation II46I I for the infiationary+coUapse 
phases. If the interior mass were plotted against the dark matter 
radius instead of the baryonic radius r;,, then the lines would be 
different, though similar. The parameters of the models and of the 
fits are listed in Table |2] 



The boundary conditions of the self-similar solutions 
are set outside the outer horizon, at a sonic point where the 
infalling baryonic plasma accelerates from subsonic to su- 
personic. Setting boundary conditions outside the horizon 
is natural, since information can propagate only inwards 
inside a black hole. The sonic point is assumed to be 
regular, meaning that the acceleration through the sonic 
point is finite and differentiable, which sets two boundary 
conditions. The accretion in real black holes is likely to 
be more complicated, but this assumption is the simplest 
physically reasonable one. 

Given the assumption of a regular sonic point, the self- 
similar solutions are characterized by four dimensionless 
free parameters: (1) the charge-to-mass ratio Q./M, of 
the black hole; (2) the mass accretion rate M,; (3) the 
equation of state w of the baryons; and (4) the dark- 
matter-to-baryon ratio pd/ Pb- Being self-similar, the so- 
lutions scale with black hole mass A/,, which increases 
linearly with time. [4^[3l also consider models where the 
baryons have a finite conductivity, and the dark matter 
has a finite interaction cross-section with the baryons, but 
the present paper assumes that such dissipative processes 
are not present. 

[49, 1(3] use a parameter 77^ in place of the dimensionless 
mass accretion rate M,; the relation between the two is 
exhibited in Appendix [Appendix E[ 
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Figure 7: Center-of-mass energy density p (thin black lines) and 
Weyl curvature scalar \C\ (purple lines almost coincident with the 
thin black lines) for the models shown in Figure \E\ labeled with 
their mass accretion rates M, = 0.003, 0.01, 0.03. The center-of- 
mass energy density and Weyl curvature inflate like the interior 
mass M. Also shown are the individual proper densities pi, of 
baryons, of dark matter, and pe of electromagnetic energy for 
the model with M, = 0.01 (to avoid confusion, only this case is 
plotted; the other models are similar). During inflation, almost 
all the center-of-mass energy p is in the streaming energy: the 
proper densities of individual components remain small. Dashed 
vertical lines mark where the outer and inner horizons would be 
in the Reissner-Nordstrom geometry with the same charge-to-mass 
Qm/M, = 0.8; mass inflation destroys the inner horizon. The density 
scale is in Planck units for a black hole of mass M. = 4 X 10'^ M0. 
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6.1. Dependence on accretion rate 

An important and surprising prediction of the analytic 
models of §21 and [5] is that the smaller the accretion rate, 
the larger the inflationary growth rate. Specifically, the in- 
flationary growth rate d\n (3 /dlnra [recall that the interior 
mass M is related to (3 by eq. ([7])] of stream a is inversely 
proportional to the "accretion rate" fib of the stream b, 
equation (I42p . where a is either of the ingoing or outgoing 
streams, and b is the other stream. This is a ferocious 
kind of instability, where the tiniest effect provokes the 
most violent reaction. 

Figure IHl shows the interior mass M as a function of 
radius for self-similar models with three different mass 
accretion rates, M, = 0.003, 0.01, and 0.03, demonstrating 
that indeed the smaller the accretion rate, the faster the 
interior mass inflates. The model with M, = 0.01 is the 
same as that illustrated in Figure S) The parameters of 
the models are listed in Table [21 

Superposed on each model in Figure [6l is the analytic 
prediction from equation (|40p for the Reissner-Nordstrom-|-infl, 
phase, and from equation (|46p for the inflation-l-collapse 
phase. The good fits confirm that the approximate ana- 
lytic models developed in §g] and [5l are a reliable guide 
to inflationary behavior, even though the models do not 



satisfy all the assumptions of the analytic approximations: 
the baryons are charged, not neutral, and the equation of 
state of the baryons is relativistic, not pressureless. 

6.2. The streams remain cold 

Figure [71 plots the energy density and Weyl curvature 
for the same set of three models as Figure [H Figure [3 
shows that the Weyl curvature scalar C and the center-of- 
mass energy density p inflate like the interior mass M . 

Figure [3 also shows the individual proper densities pb 
of baryons, pd of dark matter, and pe of electromagnetic 
energy for the model with M, = 0.01 (only this case 
is shown, to avoid confusion). The individual densities 
remain modest in spite of the fact that the center-of-mass 
energy density is inflating exponentially. This shows that 
almost all of the center-of-mass energy density is in the 
streaming energy. As in a particle accelerator, the streams 
remain cold as they are accelerated through each other. 

The fact that the proper densities of individual streams 
remain modest during inflation is consistent with the con- 
clusion of fs] that volume elements remain only mildly 
distorted during inflation, even though the Weyl curvature 
exponentiates hugely. The physical reason for the mild 
distortion is that only a tiny proper time goes by during 
inflation, as discussed in ^JUbelow. Although the tidal force 
becomes exponentially huge, it does not have enough time 
to do damage. 

7. Single fluid models 

So far, all the models considered in this paper have in- 
volved two separate fluids that stream through each other. 
However, the great majority of studies of mass inflation 
inside spherical, charged black holes have considered not 
two separate fluids, but rather the case of a massless scalar 
fiel(l usually neutral SJi^jH, M, M, HH, i, E IH 

but sometimes charged 

A massless scalar field supports waves moving at the 
speed of light, so it is plausible that mass infiation in a 
scalar field arises physically from counter-streaming be- 
tween ingoing and outgoing waves. Indeed, this was Q's 
original idea: they recognized that relativistic counter- 
streaming between ingoing and outgoing streams would 
drive the instability that they termed mass inflation, and 
they proposed that the outgoing stream would be pro- 
duced by a Price tail of outgoing gravitational radiation 
generated by gravitational collapse. A massless scalar field 
is intended to model this radiation, with the simplifying 
advantage that a massless scalar field supports spherically 
symmetric waves. 

As expounded by [s^l , the equations governing a mass- 
%s*s scalar field are the same as those of a perfect fluid 
with an ultra-hard equation of state, p — p, subject to 
the classical condition that the gradient (momentum) of 
the field is everywhere timelike. In fact the gradient of 
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Figure 8: Proper velocity V (top) of the similarity frame relative to 
the fluid, and interior mass M (bottom), as a function of radius r in 
charged single fluid models with various equations of state w = p/p. 
Curves are labeled with w. Mass inflation begins but then stalls 
in models with sound speed near but less than the speed of light, 
«; < 1. Only when the sound speed equals the speed of light, ui = 1, 
does mass inflation continue. Velocity lines are dashed where V is 
negative. Short horizontal bars mark the positions of outer and inner 
horizons, where \V\ = 1. The black hole charge-to-mass M,/Q, = 
0.8 and accretion rate M, = 0.01 are the same in all cases. 



a scalar field is not necessarily timelike everywhere, so a 
scalar field is not completely equivalent to an ultra-hard 
perfect fluid. However, an ultra-hard fluid shares with a 
scalar field the property that it supports waves moving at 
the speed of light. In what follows we suppose that the 
behavior of an ultra-hard fluid may be an adequate guide 
to the behavior of a massless scalar field. 

The models considered in this section are all self-similar, 
belonging to the same suite of self-similar models consid- 
ered in ^ but now for the case of a single fluid. 



7.1. Varying equation of state 

Is it true that mass inflation in a massless scalar field 
arises from relativistic counter-streaming between ingoing 
and outgoing waves? Support for this idea comes from 
looking at what happens in perfect fluids with various 
equations of state w — p/p. The speed of sound in a perfect 
fluid is ^/w, which equals the speed of light li w — 1, but 
is less than the speed of light if w < 1 . 

The top panel of Figure [5] shows the proper velocity V 
of the similarity frame relative to the infalling fluid for self- 
similar models with various equations of state w. These 
are single-fluid models, unlike what has been considered 
in previous sections. The fluid is charged, because in self- 
similar models the charge of the black hole is produced 
self-consistently from the charge of the accreted fluid. If 
the black hole is to be charged and hence to have an 
inner horizon where mass inflation can take place, then 
the accreted fluid must be charged. 

The velocity V in the top panel of Figure [S] equals [i^ 
the ratio C 1^*^ of the radial and time components of the 
homothetic 4-vector ^"^ (the conformal Killing vector) in 
the frame of the fiuid. The velocity V essentially defines 
the horizon structure of self-similar black holes (the fact 
that horizons can be defined unambiguously is associated 
with the presence of conformal time translation symmetry, 
that is, of self-similarity). The velocity V is sublumi- 
nal outside horizons, is equal to ±1 at horizons, and is 
superluminal between horizons. The sign of the velocity 

V provides an alternative definition of what is meant by 
ingoing and outgoing: positive V is ingoing, while negative 

V is outgoing. This definition differs from the definition of 
ingoing and outgoing adopted in the present paper, which 
is based on the sign of the radial gradient /J^, but the two 
definitions agree in the important region near the inner 
horizon. The definition adopted in the present paper is to 
be preferred because it remains valid in arbitrary spherical 
spacetimes, whereas the definition based on the sign of the 
velocity V works only if there is (conformal) time trans- 
lation symmetry. Figure [8] shows that the superluminal 
velocity V goes through infinity from positive to negative 
some way inside the horizon. This signals that the rest 
frame of the fluid changes from ingoing to outgoing, which 
happens because the charged fluid is repelled by the charge 
of the black hole. 

Figure |S] shows that mass inflation begins in models 
with sound speed near but less than the speed of light 
{w near but less than 1), but then stalls. Once inflation 
has stalled, then the fluid drops through an outgoing inner 
horizon, where the velocity is = — 1 (in the models with 
w > 0.98 in Fig. [H the numerics overflow before the inner 
horizon is actually reached). Of course, mass inflation 
would occur in these models near the inner horizon even 
if only the tiniest amount of ingoing stream were present; 
but by assumption there is only a single fluid here. 

As discussed in ^ if ingoing and outgoing streams are 
simultaneously present at an inner horizon, then necessar- 
ily the ingoing and outgoing streams must be streaming 
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Figure 9: Fluid energy density (thin blue lines) and Weyl 
curvature scalar \C\ (purple lines almost coincident with the thin 
blue lines) inside a black hole accreting a perfect fluid with an ultra- 
hard equation of state, w = p^/ = 1. Lines are labeled with their 
mass accretion rates M, = 0.001, 0.003, and 0.01. Also shown is the 
density pe of electromagnetic energy, which is the same for all the 
models. The graph illustrates that the smaller the accretion rate, 
the faster the density and curvature inflate. Dashed vertical lines 
mark where the outer and inner horizons would be in the Reissner- 
Nordstrom geometry with the same charge-to-mass Qm/M, = 0.8. 
The parameters of the models are listed in Table [2] The density 
scale is in Planck units for a black hole of mass M, = 4 X 10® Mq . 



through each other at the speed of hght. But if there is 
only a single fluid, and the speed of sound in the fluid is 
less than the speed of light, then that fluid cannot support 
both ingoing and outgoing waves at the inner horizon. If 
all streams are only ingoing, or only outgoing, then the 
mechanism of mass inflation is thwarted, and there is no 
obstacle to the streams dropping through an inner horizon, 
as seen in Figure [S) 

Figure |5] shows that only in the case where the speed 
of sound equals the speed of light, w = 1, does mass 
inflation persist. The behavior illustrated by Figure [8] is 
thus consistent with the physical idea that mass inflation 
in a massless scalar field is driven by relativistic counter- 
streaming between ingoing and outgoing waves. 

7.2. Dependence on accretion rate 

Figure M shows the energy density and Weyl curva- 
ture scalar |C| for self-similar models accreting an ultra- 
hard fluid, w = p^/ — I, at three different mass accre- 
tion rates M, = 0.001, 0.003, and 0.01. The model with 
M, = 0.01 was already featured in Figure [8] 

Figure [5] shows that the inflationary growth rate is 
faster for smaller accretion rates, in agreement with ex- 
pectations from the analytic models of §21 ^nd [SI As 
illustrated in Figure IIOI the inflationary growth rate is 
inversely proportional to the mass accretion rate, except 
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Figure 10: Inflationary growth rate |dln/3/dlnrj as a function of 
mass accretion rate M, in self-similar models in which an ultra-hard 
fluid accretes on to a black hole with charge-to-mass Q./Af, = 0.8. 
The inflationary growth rate is inversely proportional to the mass 
accretion rate, equation II50II . except at the largest accretion rates. 
The curve terminates at a maximum accretion rate of M, = 0.0784. 



at the largest accretion rates. For accretion rates away 
from the maximum, the curve shown in Figure [10] fits to 

^ . . (50) 

d\nr M, 

The maximum possible accretion rate is constrained by 
the requirement of appropriate boundary conditions at the 
sonic point ^49j . For an ultra-hard fluid accreting on to a 
black hole of charge-to-mass Q,/M, = 0.8, the maximum 
accretion rate is M, = 0.0784. 

Figure [HI shows that (at least for constant accretion 
rate, as is necessarily true in self-similar models) the ultra- 
hard fluid shows a single power-law inflationary phase, 
as opposed to two distinct phases of inflation followed 
by collapse. The ultra-hard fluid does collapse to zero 
radius (if quantum gravity does not intervene) , but not as 
a separate phase. The interior mass tends to infinity at 
zero radius, as opposed to reaching a finite value. 

A prominent difference between the ultra-hard fluid 
and the two-stream models considered in previous sections 
is that the ultra-hard fluid always maintains a large trans- 
verse pressure. Of course this is built in to the equation 
of state p = p oi the ultra-hard fluid. But it does contrast 
with the two-stream models, where inflation accelerates 
ingoing and outgoing streams in opposite directions, which 
produces an exponentially growing radial pressure in the 
center-of-mass frame, but no transverse pressure to go 
with it. Which is a correct, a large or a small transverse 
pressure? The authors' inclination is to suspect that the 
two-stream model, with small transverse pressure, is likely 
to be a closer model of reality, on the grounds that the feed- 
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back loop that drives inflation involves only radial accel- 
erations and radial components of the energy-momentum 
tensor, with no dependence on transverse pressure. Absent 
some physical mechanism to maintain isotropy, it is hard 
to see how the transverse pressure could keep up with the 
exponentially growing radial pressure. 

8. The far future? 

The Penrose diagram of a Reissner-Nordstrom or Kerr- 
Newman black hole indicates that an observer who passes 
through the outgoing inner horizon sees the entire future 
of the outside universe go by. In a sense, this is "why" the 
outside universe appears infinitely blueshifted. 

This raises the question of whether what happens at 
the outgoing inner horizon of a real black hole indeed de- 
pends on what happens in the far future. If it did, then the 
conclusions of previous sections, which are based in part 
on the proposition that the accretion rate is approximately 
constant, would be suspect. A lot can happen in the far 
future, such as black hole mergers, the Universe ending 
in a big crunch. Hawking evaporation, or something else 
beyond our current ken. 

Ingoing and outgoing observers both see each other 
highly blueshifted near the inner horizon. An outgoing 
observer sees ingoing observers accreted to the future of 
the time that the outgoing observer fell in, while an ingo- 
ing observer sees outgoing observers accreted to the past 
of when the ingoing observer fell in. If the streaming 
4- velocity between ingoing and outgoing streams is w^, 



equation (jC.2|) . then the proper time dri, that elapses on 
stream b observed by stream a equals the blueshift factor 
u^jj times the proper time dTa experienced by stream a. 



dn = ul^ dTa 



(51) 



8.1. Inflationary phase 

A physically relevant timescale for the observing stream 
a is how long it takes for the blueshift to increase by one e- 
fold, which is dTa/dlnul^. During the inflationary phase, 
stream a sees the amount of time elapsed on stream b 
through one e-fold of blueshift to be 



dTa 



d\nul^ 



(52) 



The right hand side of equation ([5^ is derived from u^^ « 
2\uaUb\, \dra/dTa\ = \l3a,t\ ~ Pua, Ta ~ r_ , and the ap- 
proximations (|^T|) and (|42l) valid during the inflationary 
phase. The constants Cf, and A on the right hand side 
of equation (|52p are typically of order unity, while r_ 
is the radius of the inner horizon where mass inflation 
takes place. Thus the right hand side of equation (1521) 
is of the order of one black hole crossing time. In other 
words, stream a sees approximately one black hole crossing 
time elapse on stream b for each e-fold of blueshift. The 
estimate ([5^ is compatible with formula (1) of f57^, who 



considered an outgoing infaller irradiated by cosmic mi- 
crowave background photons, truncated when the blueshifted 
photons reached Planck energy. 

For astronomically realistic black holes, exponentiat- 
ing the Weyl curvature up to the Planck scale will take 
typically a few hundred e- folds of blueshift, as illustrated 
for example in Figure [T] Thus what happens at the inner 
horizon of a realistic black hole before quantum gravity 
intervenes depends only on the immediate past and future 
of the black hole - a few hundred black hole crossing times 
- not on the distant future or past. This conclusion holds 
even if the accretion rate of one of the ingoing or outgoing 
streams is tiny compared to the other, as considered in 



From a stream's own point of view, the entire inflation- 
ary episode goes by in a flash. At the onset of inflation, 
where /3 goes through its minimum, at /iaW^ — iJ^bul — A 
according to equation ([57]) . the blueshift is already 
large 

2A 



rmr ^^^^ 

thanks to the small accretion rates /Za and ^b- During the 
first e-fold of blueshift, each stream experiences a proper 
time of order y/JI^JIb times the black hole crossing time, 
which is tiny. Subsequent e-folds of blueshift race by in 
proportionately shorter proper times. The time is so short 
that volume elements of each stream are little distorted 
during inflation despite the huge tidal force, as first pointed 
out by [|. 

8.2. Collapse phase 

The time to reach the collapse phase is another matter. 
According to the estimates in ^ ^5.'6\ and 15.41 reaching the 
collapse phase takes of order ^ 1/ fXa e-folds of blueshift, 
where jjLa is the larger of the accretion rates of the ingoing 
and outgoing streams, equation (|49)) . Thus in reaching the 
collapse phase, each stream has seen approximately l//^a 
black hole crossing times elapse on the other stream. But 
l/jJLa black hole crossing times is just the accretion time 
- essentially, how long the black hole has existed. This 
timescale, the age of the black hole, is not infinite, but it 
can hardly be expected that the accretion rate of a black 
hole would be constant over its lifetime. 

If the accretion rate were in fact constant, and if quan- 
tum gravity did not intervene and the streams remained 
non-interacting, then indeed streams inside the black hole 
would reach the collapse phase, whereupon they would 
plunge to a spacelike singularity at zero radius. For ex- 
ample, in the self-similar models illustrated in Figures H] 
and [3 outgoing baryons hitting the central singularity see 
ingoing dark matter accreted a factor of two into the future 
(specifically, for ingoing baryons and outgoing dark matter 
hitting the singularity, the radius of the outer horizon 
when the dark matter is accreted is twice that when the 
baryons were accreted; the numbers are 2.11 for Af, = 
0.03, 1.97 for M. = 0.01, and unknown for M. = 0.003 
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because in that case the numerics overflow before the cen- 
tral singularity is reached). The same conclusion applies 
to the ultra-hard fluid models illustrated in Figure [SI if 
the accretion rate is constant, then outgoing streams see 
only a factor of order unity or a few into the future before 
hitting the central singularity. 

If on the other hard the accretion rate decreases suf- 
ficiently rapidly with time, then it is possible that an 
outgoing stream never reaches the collapse phase, because 
the number of e-folds to reach the collapse phase just keeps 
increasing as the accretion rate decreases. By contrast, an 
ingoing stream is always liable to reach the collapse phase, 
if quantum gravity does not intervene, because an ingoing 
stream sees the outgoing stream from the past, when the 
accretion rate was liable to have been larger. 

However, as already remarked in in astronom- 

ically realistic black holes, it is only for large accretion 
rates, such as may occur when the black hole first collapses 
(M. > 0.01 for the models illustrated in Fig. [T]), that the 
collapse phase has a chance of being reached before the 
Weyl curvature exceeds the Planck scale. 

To summarize, it is only streams accreted during the 
first few hundred or so black hole crossing times following 
a black hole's formation, or following an exceptionally high 
accretion event such as a black hole merger, that have a 
chance of hitting a central spacelike singularity. Streams 
accreted at later times, whether ingoing or outgoing, are 
likely to meet their fate in the inflationary zone at the 
inner horizon, where the Weyl curvature exponentiates to 
the Planck scale and beyond. 

8.3. Null singularity on the Cauchy horizon? 

[4^ has proved a number of mathematical theorems 
that establish that a null singularity forms on the Cauchy 
horizon of a charged spherical black hole accreting a mass- 
less scalar field. The situation envisaged by the theorems 
is that of a black hole that collapses and thereafter remains 
isolated. The collapse generates an outgoing Price tail of 
radiation. The theorems assume that the outgoing Price 
radiation falls off sufficiently rapidly along outgoing null 
geodesies, and [2^ have proved that the required condition 
on the Price radiation holds for an isolated spherical black 
hole accreting a massless scalar field. The theorems con- 
firm the several analytic and numerical studies that have 
found a null singularity on the Cauchy horizon 

The conclusion is completely consistent with the results 
of the present paper. However, as regards real astronomi- 
cal black holes, the question of whether a null singularity 
forms is academic. As emphasized by Q and discussed in 
^8.1[ quantum gravity is likely to intervene. If quantum 
gravity is set aside, then, as discussed in ^18. 21 the appear- 
ance of a null singularity depends on events happening in 
the indefinite future. 

M, E] find erically that a null singularity forms 



sufficiently rapidly, the required degree of rapidity depend- 
ing on the parameters of the problem, such as the charge- 
to-mass ratio of the black hole. If too much scalar field con- 
tinues to be accreted, then no null singularity forms, and 
the field collapses to a central singularity. These results are 
consistent with the arguments of the present paper. Again 
however, we emphasize that the appearance of a null singu- 
larity depends on events happening in the indefinite future, 
and neglects quantum effects that undoubtedly will be 
important, whether those arise from super-Planckian den- 
sities and curvatures, or from super-Planckian collisions 
between ingoing and outgoing streams, or from Hawking 
evaporation of the black hole, or from pair creation near 
the inner horizon. 



58| have shown that in the simplified case of a 1-1-1- 



dimensional charged black hole, if the effects of pair cre- 
ation of charged particles are taken into account, then the 
result is collapse to a spacelike singularity rather than 
a null singularity on the Cauchy horizon. The result is 
consistent with the argument of the present paper that 
as long as there is any source that continues to replenish 
ingoing and outgoing streams near the inner horizon, the 
ultimate result will be collapse to a spacelike singularity. 
The results of ^58] suggest that even without any direct ac- 
cretion, pair creation provides a sufficient source of ingoing 
and outgoing streams. 

9. Collision rate 

It has been assumed throughout this paper that ingoing 
and outgoing streams stream through each other without 
interacting. The aim of this section is to check the validity 
of that assumption. The conclusion is that the assumption 
is good at typical astronomically low accretion rates, but is 
likely to break down at higher accretion rates. Moreover 
the assumption may well break down at center-of-mass 
collision energies exceeding the Planck mass. We hope to 
explore the consequences of interaction between streams 
in a subsequent paper. 

For definiteness, the models considered in this section 
are the baryon-plus-dark-matter models illustrated in Fig- 
ure[71 As a practical matter, it is unlikely that, as assumed 
by the models, the ingoing stream would consist only of 
dark matter and the outgoing stream only of baryons. 
More realistically, both streams would contain baryons, 
and baryon-baryon collisions would dominate baryon-dark 
matter collisions, at least at low energies. However, it is 
convenient to refer to the two streams as "baryon" and 
"dark matter" streams. 

Each baryon in the black hole particle accelerator sees 
a fiux ridu of dark matter particles per unit area per unit 
time, where rid = Pd/'nid is the proper number density 
of dark matter particles in their own frame, and u = 
uj^j, is the radial component of the proper streaming 4- 
velocity, the 7U, of the dark matter through the baryons 



only if the scalar field set up outside the horizon falls off brevity, in this subsection the streaming velocity 



equation (|C.2I) . is written simply as it; this is not the same 
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Figure 11: Collision rate of the black hole particle accelerator per 
e-fold of velocity u (meaning 711), expressed in units of the inverse 
black hole accretion time M,/M,. The models illustrated are the 
same as those in Figure [7] The curves are labeled with their mass 
accretion rates: M, = 0.03, 0.01, 0.003, and lO"!^. Stars mark 
where the center-of-mass energy of colliding baryons and dark matter 
particles exceeds the Planck energy, while disks show where the Weyl 
curvature scalar C exceeds the Planck scale. 



as the u defined by equation ©]. The 7 factor in u is 
the relativistic beaming factor: all frequencies, including 
the collision frequency, are speeded up by the relativistic 
beaming factor 7. As the baryons accelerate through the 
collider, they spend a proper time interval dr/dlnu in 
each e-fold of velocity u. The number of collisions per 
baryon per e-fold of u is the dark matter flux (pd/mci)u, 
multiplied by the time dr/dlnu, multiplied by the collision 
cross-section a. 



number of collisions 
baryon x e-fold of u 



Pd dr 
rud d In u 



(54) 



During the inflationary phase, the collision rate coefficient 
PdU dr/dlnu that goes into equation (|54p is, from equa- 
tion ([5^ and the definition ([M)) of p,d, 



PdU 



dr 
dlnu 



87rAr_ ^ M. ' 



(55) 



in which the proportionality on the rightmost side follows 
from fid oc M, and r_ cx Af,. Figure [TT] shows, for 
several different mass accretion rates M, , the collision rate 
PdU dr/dlnu of the black hole collider, expressed in units 
of the inverse black hole accretion time M,/M,. In the 
units c = G = 1 being used here, the mass of a baryon 
(proton) is 1 GcV w 10^^^ m. If the cross-section cr is 
expressed in units of femtobarns (Ifb = 10"'*'^ m^), which 
is about a weak interaction cross-section, then the number 
of collisions ([M)) per baryon per e-fold of collision velocity 



number of collisions 
baryon x e-fold of u 

( M,/M, 
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(1010 yr)- 



' PdU dr / dlnu 
, 0.03 A/,/M. 



(56) 



The number (|56p of collisions per baryon is small. This is 
like a particle accelerator, where collisions are rare. If the 
dark matter stream is dominated by baryons, as may well 
be true in realistic models, then the collision rate ([551) rises 
to of order unity for low-momentum-transfer electromag- 
netic or strong collisions, where cross-sections approach 
one barn. However, a few low-momentum-transfer colli- 
sions will not spoil the conclusion that to a good approxi- 
mation the counter-streaming streams are non-interacting, 
at least in the realm of "known physics," at sub-Planckian 
center-of-mass collision energies. 

The collision rate (l55|) is scaled to the typical accretion 
timescale M,/Af, ^ lO^^yr of astronomical black holes. 
The collision rate would be higher during episodes of high 
accretion. At sufficiently high accretion rates, collisions 
between streams would affect the streams, and the as- 
sumption of non-interacting streams would break down. 

Stars in Figure [TT] mark where the center-of-mass en- 
ergy y/mbiTidU of colliding baryons and dark matter parti- 
cles, taken to have masses nib = 1 GeV and rud = 300 GeV, 
hits the Planck energy. This occurs well before the cur- 
vature hits the Planck scale. If cross-sections increase 
sufficiently rapidly at super-Planckian center-of-mass ener- 
gies, then again collisions between streams would become 
important, and the assumption of non- interacting streams 
would fail. 

The total number of collisions taking place in the black 
hole per unit (external) time equals the rate ([55]) per baryon 
multiplied by rate at which the black hole is accreting 
baryons, which is approximately equal to the mass accre- 
tion rate M, divided by the mass rrib per baryon: 



number of collisions 
time X e-fold of u 



M, Pd dr 

^^T^ • ^'^) 

TOfc nid dlnu 



Numerically, the number of collisions (|57p per unit (exter- 
nal) time is 



number of collisions 



time X e-fold of u 

M, \ I M,/M, 



a \ /300GcV^ 




Ifb/ \ Tfibrrid 
Pdudr / dlnu^ 



(10iOyr)"M V0.03M./A/. 



(58) 



The collision rate ([58]) is scaled to the accretion rate of 
the Milky Way black hole, SJ2] Current-era particle ac- 
celerators are proud of delivering inverse femtobarns of 
luminosity during their lives. Equation (|58p shows that the 
accelerator inside the Milky Way black hole is delivering 
of order 10'^^ inverse femtobarns per year in each e-fold of 
collision velocity up to the Planck energy and beyond. 
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10. Summary 

The purpose of this paper has been to give a clear ac- 
count of the physical causes underlying the mass inflation 
instability first proposed by Q . The arguments have been 
restricted to the case of spherical, charged black holes, 
which has been the case considered by most studies of mass 
inflation to date. It is generally thought that charge may 
be a satisfactory surrogate for spin, since charged black 
holes have inner horizons like rotating black holes. 

The elements of mass inflation can be seen already in 
the two Einstein equations valid for arbitrary spheri- 
cally symmetric spacetimes. The first (15a)) of these equa- 
tions relates the acceleration to the gravitational force, as 
measured by a freely-falling observer. Acceleration here 
means the proper rate of change of the proper radial veloc- 
ity /3t = dtr measured by the observer. The equation (I5a|) 
differs from the corresponding Newtonian equation in that 
the gravitational force is sourced not only by the interior 
mass, but also by radial pressure. A second essential 
difference compared to Newtonian gravity is that gravity in 
spherical spacetimes is governed by two equations, not one. 
The radial velocity /3t is one component of a 4-vector, the 
radial 4-gradient (3m = = {/3t, /3r, 0, 0}. The Einstein 
equation ([8b)) for the second component Pr is also central 
to producing mass inflation. 

As originally argued by Q, mass inflation is produced 
by relativistic counter-streaming between ingoing and out- 
going streams just above the inner horizon. The present 
paper has shown that the counter-streaming is driven by 
the gravitational force produced by the streaming pressure 
and energy flux. This is why mass inflation is exponential: 
the streaming pressure and flux increase the gravitational 
force, which accelerates the streams faster through each 
other, which increases the streaming pressure and flux, 
which increases the gravitational force, and so on. 

A crucial ingredient of the process is that the gravi- 
tational force acts in opposite directions for ingoing and 
outgoing streams: it accelerates ingoing streams towards 
the black hole, and outgoing streams away from the black 
hole. [l6| explain this conundrum by pointing out that 
the gravitational force is always inwards, in the direction 
of smaller radius, and that the direction of smaller radius 
is towards the black hole for ingoing streams, and away 
from the black hole for outgoing streams. The Einstein 
equation (|8b| shows that this picture is essentially correct, 
although as argued in the next paragraph it is too simplis- 
tic. A stream is ingoing or outgoing depending on whether 
the proper radial gradient /^^ = drr measured in that frame 
is positive or negative, and the Einstein equation ()8bp 
shows that increases or decreases depending on whether 
the energy flux in that frame is positive or negative. The 
result is that gravity tends to drive a positive (ingoing) /Sr 
to become more positive, and a negative (outgoing) /Sr to 
become more negative. 

However, this picture that the gravitational force points 
in the direction of smaller radius, which is in opposite di- 



rections for ingoing and outgoing streams, does not explain 
why mass inflation can come to an end, as found in the 
simple analytic models constructed in §2laiid[S] The more 
complete picture painted by the Einstein equations ([5]) is 
that the 4-vector nature of /?,„, whose time component is 
the radial velocity l3t, is central. The ingoing and outgoing 
streams are accelerated through each other only if the 
change in the radial component f3r exceeds the change 
in the time component /3t, as illustrated in the middle 
panel of Figure [1] If the change in [3r is less than that 
in (3t, then the streaming velocity decreases, as illustrated 
in the bottom panel of Figure [TJ If this happens, then 
mass inflation comes to an end. 

Sections|3]and[5]build on the qualitative picture of 33|to 
develop simple approximate analytic models that capture 
quantitatively how inflation is ignited, exponentiates, and 
then ends. The flrst of the two sections, 21 considers 
the simple but unrealistic example of symmetrically equal 
ingoing and outgoing streams, which 4^ 47,|33i] has pre- 
viously called the homogeneous approximation (see 
The second of the two sections, ^ considers the realistic 
case of unequal ingoing and outgoing streams. The reader 
interested in understanding mathematically how Einstein's 
equations ([5]) imply mass inflation is encouraged to try the 
problem in Appendix jAppendix B| 

The predictions of the analytic approximations are borne 
out by comparison to self-consistent self-similar solutions. 
Section[6|considers self-similar solutions in which the black 
hole accretes two separate streams, "baryons" and "dark 
matter," while 53 considers self-similar solutions in which 
the black hole accretes just a single fluid. The single-fluid 
models include the case of a perfect fluid with an ultra- 
hard equation of state p = p, which can be taken as a 
model of a massless scalar field, subject to the constraint 
that the momentum of the field is everywhere timelike (Hoj . 

In the single-fluid models considered in ^Zl it is shown 
that mass inflation occurs only in the case of the ultra- 
hard fluid, whose sound speed equals the speed of light. If 
the equation of state w ~ p/ p is near but not equal to 1, 
so that the sound speed ^Jw is near but not equal to the 
speed of light, then mass inflation begins, but then stalls, 
whereupon the outgoing charged fluid drops through the 
outgoing inner horizon. The behavior is consistent with 
the notion that mass inflation in a fluid that supports 
waves moving at the speed of light, such as a massless 
scalar field, results from relativistic counter-streaming be- 
tween ingoing and outgoing waves, as originally suggested 

by i- 

Two possibly surprising predictions of the analytic so- 
lutions are: first, that the smaller the accretion rate, the 
faster infiation grows; and second, that infiation eventually 
comes to an end, whereupon the ingoing and outgoing 
streams collapse to a spacelike singularity at zero radius. 
Both predictions are confirmed in the self-consistent solu- 
tions, in both two-fiuid and single-fluid cases. 

The fact that a smaller accretion rate causes faster 
inflation means that the inflationary instability is difficult 
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to avoid. How is the beast to be tamed if the tiniest effect 
provokes the most violent reaction? 

One possible way to avoid inflation is to introduce a 
large amount of dissipation j59ij , fast enough to transport 
charge and angular momentum so that the black hole be- 
comes neutral and non-spinning towards its center, elimi- 
nating the inner horizon where inflation is ignited. Given 
that angular momentum transport is a rather weak process 
(60| , we are inclined to suspect that real rotating black 
holes do not dissipate all their spin, and that inflation 
does occur in reality. 

The results of this paper are consistent with all the 
previous literature on the mass inflation instability. It 
might seem that there is a conflict in that we do not 
confirm the often-stated conclusion that the generic conse- 
quence of inflation is a weak null singularity on the Cauchy 
horizon. The apparent conflict arises because most studies 
have considered the situation of a black hole that collapses 
and thereafter remains isolated, whereas we have taken 
the point of view that a real black hole is never isolated. 
As long as a black hole continues to accrete, and thereby 
to generate ingoing and outgoing streams near the inner 
horizon, the inflationary growth rate remains finite, and 
there is no null singularity. Even in the absence of direct 
accretion, quantum-mechanical pair creation provides a 
source of ingoing and outgoing streams (Hsj . preventing 
a null singularity. 

For typical astronomically low accretion rates, inflation 
drives the Weyl curvature and the center-of-mass energy 
density and pressure far above the Planck scale before sig- 
nificant tidal distortion occurs, as first pointed out by [8|. 
Physically, even though tidal forces grow super-Planckian, 
the timescale for those forces to act is so short that volume 
elements remain little distorted. Consequently, as already 
pointed out by Q , quantum gravity in some form is likely 
to intervene before any classical general relativistic singu- 
larity is reached. 

If on the other hand the accretion rate is large enough, 
as happens when the black hole first collapses, and during 
events of high accretion such as a black hole merger, then 
the streams collapse to a central singularity before infla- 
tion drives the curvature above the Planck scale. Such 
a crushing singularity is subject to non-radial Belinski 
Khalatnikov-Lifshitz instabilities 



not been explored in this paper; see 




55| . which have 
for a review. 



11. Final remarks 

11.1. Misconceptions 

As a pedagogical device, it is helpful to confront a num- 
ber of misconceptions about what happens inside black 
holes. The following remarks are predicated on the as- 
sumption that what happens in a spherical, charged black 
hole, as considered in this paper, is a reliable guide to what 
happens in a realistic rotating black hole. 

Misconception 1. A popular story is that matter that 
falls inside a black hole falls to a central singularity, a point 



of infinite curvature, where space and time stop. This mis- 
conception is based on the prototype of the Schwarzschild 
geometry, bolstered by a simplistic interpretation of the 
Penrose-Hawking 6l|, |62| singularity theorems. In reality, 
as long as a black hole has any spin, as is certain to be true 
(the converse is a set of measure zero), and except during 
epochs of high accretion, such as within a few hundred 
black hole crossing times of the initial collapse or of a 
black hole merger, then matter falling into the black hole 
will reach its nemesis in the inflationary region at the 
inner horizon, which will be at a macroscopic radius (some 
fraction of the radius of the outer horizon) as long as 
the spin of the black hole is appreciably different from 
zero. Except during epochs of high accretion, inflation 
exponentiates the curvature and proper density beyond 
the Planck scale before collapse to a central singularity oc- 
curs. What happens after the curvature and density have 
reached Planck scale depends on quantum gravitational 
processes that have yet to be explored in the literature. 

Misconception 2 (^S]). What happens inside a black 
hole depends on what happens in the infinite future. This 
misconception comes from looking at the Penrose diagram 
of the Reissner-Nordstrom or Kerr-Newman geometry. In 
reality, during inflation, an outgoing observer sees approxi- 
mately one black hole crossing time elapse on ingoing mat- 
ter for each e-fold of blueshift (an outgoing observer sees 
ingoing matter accreted from the future, while an ingoing 
observer sees outgoing matter accreted in the past). After 
a few hundred e-folds, the curvature will have exponen- 
tiated to the Planck scale, whereupon quantum gravity 
presumably intervenes. Far from depending on the infinite 
future, what happens to an outgoing or ingoing stream 
undergoing inflation depends on events only a few hundred 
black hole crossing times into the future or past. 

Misconception 3 ( j)8.3p . Mass inflation produces a null 
singularity on the Cauchy horizon. As proven by 42| , 
confirming earlier studies H, E, El, [III S, IH, [ijra 
null singularity occurs on the Cauchy horizon (outgoing 
inner horizon) of a black hole that collapses and thereafter 
remains isolated for ever. However, a real black hole is 
never isolated, and for ever is a long time. It is true 
that an (outgoing or ingoing) observer who falls into an 
astronomical black hole is likely to meet their fate at the 
inner horizon where inflation exponentiates the curvature 
to the Planck scale and beyond. But that fate is something 
to do with quantum gravity, not a general relativistic null 
singularity. 



11.2. The black hole particle accelerator 

We conclude this paper on a note of wonder. If the 
story told in this paper is true, then Nature has devised a 
most remarkable accelerator of extraordinary power inside 
black holes. The accelerator is powered by gravity, and 
it feeds on itself: the higher the streaming energy, the 
stronger the gravitational acceleration. The accelerator 
is ingeniously constructed so that the smaller the initial 
streams, the more rapidly the acceleration exponentiates. 
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ensuring that nothing can escape it. What does Nature do 
with such a machine? 
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Appendix A. Spherically symmetric spacetime 

The orthonormal tetrad formalism introduces at each 
point of spacetime a tetrad, a set of orthonormal axes 7m- 
The tetrad is related to the basis of coordinate tangent 
vectors by the vierbein e™'' and its inverse e™^ 

7m = Gm^g^ , Qf, = e"\,'Jrn ■ (A.l) 

By construction, the scalar products of the orthonormal 
tetrad axes fm constitute the Minkowski metric Tymn , while 
by definition, the scalar products of the tangent basis 
constitute the coordinate metric g^i, 

The line-element (O encodes not only a metric but also an 
inverse vierbein e™^, through 

e*^ (Ixf^ — dt/a , 



and the corresponding vierbein e^'^ is 



e'-^dx^ = {l/l3r)idr-Ptdt/a) 
e'^^dx" :^rdd , 



ef^ dxf" = r sin 61 dcf) 



Explicitly, the inverse vierbein e 
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The tetrad- frame connection coefficients Tkmn are 





= ht , 


(A.6a) 






(A.6b) 
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where j3m are defined by equation and hm are defined 
in terms of the vierbein coefficients by 



/li = — 9r In a , 

hr = -Pt^ h Otln/Sr 

or or 



(A.7a) 
(A.7b) 



A person at rest in the tetrad frame has, by definition, 
tetrad- frame 4- velocity u*^ = {1, 0, 0, 0}, and the covariant 
derivative DmU^ of the radial component of the tetrad- 
frame 4-velocity of such a person is 



DmU'' 



(A.; 



Equations (|A.8I) reveal that ht and hr represent respec- 
tively the proper radial acceleration (minus the gravita- 
tional force), and the "Hubble parameter" of the radial 
flow, experienced by a person at rest in the tetrad frame. 
Restricted to the i-r-plane, hm = {ht,hr} constitute a 
tetrad-frame 2-vector. 

The non- vanishing components of the tetrad-frame Rie- 
mann tensor Rkimn are 



Rtete 
Rtdre - 

RrffrO - 



Rtrtr 

- Rt<f>t4> 

- Rt4>r4> 
Rr(j)r(f> 
if, 



— Drht — Dthr 



(A.9a) 
(A.9b) 

(A.9c) 

(A.9d) 

(A.9e) 

The non-vanishing components of the tetrad-frame Ein- 
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stein tensor G*'™ are 



M 



r \ r'^ J 

r \ y 
1 



(A.lOa) 
(A.lOb) 
(A.lOc) 



(A.lOd) 

In spherically symmetric spacetimes, the only distinct non- 
vanishing component of the Weyl curvature tensor is the 
scalar (spin-0) component, the Weyl curvature scalar C 

"^/^ ^ 



6 



(G" - G'"' + G"") 



(A.ll) 



4-gradient = = {dtr, drr, 0, 0} in the frame of any 
radially moving orthonormal tetrad is 



D,I3, = 4»r/ 



(B.la) 
(B.lb) 



where Dt is the covariant time derivative, p is the radial 
pressure, / is the radial energy flux, and M is the interior 
mass defined by 



2M 
r 



1 = 13^ 



(B.2) 



Appendix B.l. Freely-falling stream 

Consider a stream of matter that is freely falling ra- 
dially inside the horizon of a spherically symmetric black 
hole. Let u be the radial component of the tetrad-frame 4- 
velocity u™ of the stream relative to the "no-going" frame 
where /3r = (the frame of reference that divides ingoing 
frames /3r > from outgoing frames (ir < 0): 



Thanks to the Bianchi identities, the Einstein tensor au- 
tomatically satisfies covariant conservation 



D^G-" 



, 



(A.12) 



enforcing conservation of energy-momentum. The two non- 
vanishing conservation equations are those for energy and 
radial momentum, which, written out in full, are 



dt + '^-^ + hA G" 
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dr + — + 2ht] G*'' 



/i, G'^'- + — G^^ = , (A. 13a) 



dt + '^ + 2k\ g*'- +(dr • ^^"^ 
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ht 1 G'^'' 



h, G" - ^ G« 



(A.13b) 



Since each of the two conservation equations (jA.13[) in- 
volves each of the 4 distinct non-vanishing components 
G", G*'', G'''^, and G^^ of the Einstein tensor, it follows 
that any two of the Einstein equations may be dropped in 
favor of the two equations of energy-momentum conserva- 
tion. 

Appendix B. Problem 

This problem appeared on the final exam of the first 
author's graduate course on general relativity in Spring 
2008. Some of the equations below repeat equations else- 
where in this paper, but they are left as is so that the 
problem remains self-contained. 

The mechanism of mass inflation 

Einstein's equations in a spherically symmetric space- 
time imply that the covariant rate of change of the radial 



u™ = {-/3t//3, -/3,//3, 0, 0} = Wl + u^, u, 0, 0} . (B.3) 

Note that is negative inside the horizon. The time 
component = — /3t//3 = vT+z? of the 4- velocity is 
positive (as it should be for a proper 4- velocity), while the 
radial component u = = —I3r/f3 of the 4- velocity is 
negative ingoing, positive outgoing. Show that along the 
worldline of the stream 
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[Hint: If the stream is freely falling, then the proper time 
derivative dt in the tetrad frame of the stream equals 
the covariant time derivative Df. Thus the proper rates 
of change of ln/3 and Ium with respect to Inr along the 
worldline of the stream are 



dln^ dt\nl3 



din? 



d In r 9t In r ' d In r 
These can be evaluated through 



dt Inu 
dt Inr 
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dt ln/3 = A ln/3 = T7^Dt(3' = TT^Dti^t " 
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(B.6c) 



with Einstein's equations (jB.ip substituted into equations 
([RHa)) and ([RBbl) .] 
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Appendix B.2. Equal ingoing and outgoing streams 

Consider the symmetrical case of two equal streams 
of radially ingoing (/?,- > 0) and outgoing {j3r < 0) neu- 
tral, pressureless, non- interacting matter ("dust"), each of 
proper density p in their own frames, freely- falling into a 
charged black hole. Show that 



dln/3 
dlnr 
dlnzi 
dlnr 



where 



Hence conclude that 

dln^ -\ + 13"^ + ^lu^ 



dlxi'i 



X ~ P'^ + fi + p.u'^ 



(B.7a) 
(B.7b) 

(B.8) 
(B.9) 



[Hint: The assumption that the streams are neutral, pres- 
sureless, and non-interacting is needed to make the streams 
freely- falling, so that equations (jB.4p are valid. The pres- 
sure p in the tetrad frame of each stream is the sum of the 
electromagnetic pressure Pe and the streaming pressure Ps 



p^Pe+Ps ■ 

The electromagnetic pressure Pe is 



Pe 
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(B.IO) 



(B.ll) 



with Q the charge of the black hole, which is constant 
because the infalling streams are neutral. The streaming 
pressure Ps that each stream sees is 



Ps 



r\2 



(B.12) 



where the streaming 4- velocity u™ between the two streams 
is the 4- velocity of the observed stream Lorentz-boosted by 
the 4- velocity of the observing stream (the radial velocities 
of the observed and observing streams have opposite 
signs) 



(u*)2 + (u'')2 = 1 + 2u2 



(B.13a) 
(B.13b) 



The energy flux / in the tetrad frame of each stream is the 
streaming flux fs 



p r t r 

f = fs= PU^U^ 



(B.14) 



You should find that the combinations of streaming pres- 
sure and flux that go into equations (|B.4p are 



Ps + ^fs = '^pu^ , 
Pt 

Ps + ^fs = -2p(l + 



(B.lSa) 
(B.15b) 



Appendix B.3. Reissner- Nordstrom phase 

If the accretion rate is small, then initially the stream 
density p is small, and consequently p is small. Argue that 
in this regime equation (jB.9p simplifies to 



dln^ + 



dlnw A-/32 
Hence conclude that 



(B.16) 



(B.17) 



where C is some constant set by initial conditions (gener- 
ically, C will be of order unity). 

Appendix B.4-- Transition to mass inflation 

Argue that in the Reissner-Nordstrom phase, /? be- 
comes small, and u grows large, as the streams fall to 
smaller radius r. Argue that in due course equation (jB.Qp 
becomes well-approximated by 



d\iif3 — A + /iu^ 



dlnu A + pu^ 



(B.18) 



Treating A and p as constants (which is a good approxi- 
mation) , show that the solution to equation (jB.lSp subject 
to the initial condition set by equation (IB.17P is 



/3 



C{X + pu^) 
Xu 



(B.19) 



[Hint: A is positive. In the Reissner-Nordstrom solution, 
/3 would go to zero at the inner horizon.] 

Appendix B.5. Sketch 

Sketch the solution (jB.lQp . plotting u against j3 on 
logarithmic axes. Mark the regime where mass infiation 
is occurring. 

Appendix B.6. Inflationary growth rate 

Argue that during mass infiation the infiationary growth 
rate d In /S/d In r is 



dln/3 
dlnr 



A2 



(B.20) 



Comment on how the inflationary growth rate depends on 
accretion rate (on p). 

Appendix C. Multiple streams 

Let there be multiple neutral, pressureless, radially 
freely- falling streams. Analogously to equation (jB.Sp . let 
Ua denote the radial component of the 4-velocity of 
stream a relative to the no-going frame: 



C = {-/3a,t//3,-/3a,r//3,0,0} = {v/r+^,Wa,0,0} 



(C.l) 
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The velocity Ua is negative ingoing, positive outgoing. Let 
denote the 4-velocity of stream b relative to stream a 



'^ba — '^b'^*'a 



UuU„ 



UuU„ 



(C.2) 



The proper pressure seen by stream a is the sum of the 
electromagnetic pressure Pe, equation (|B.11I) . which is the 
same for every stream, and the streaming pressure Pa, 
which is a sum over the contributions from all other streams 
b: 



Pa 



b 



Pbiula) 



The energy flux fa seen by stream a equals the streaming 
energy flux, which is again a sum over the contributions 
from all other streams b: 



PbUbaUba 



(C.4) 



Equations (|B.4|) remain valid individually for each freely- 
falling stream. The combinations of streaming pressure 
and flux that go into equations (|B.4|) are 
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fa = ^pbiuiy 
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(C.5a) 
(C.5b) 



In the particular case of two streams, one ingoing and 
one outgoing, that are streaming relativistically through 
each other, as is the relevant case during mass inflation. 



equations (|C.5|) simplify to 



Pa 



Pa 



Pa,t 



fa « 2pbul , 

fa W -2pbul 



(C.6a) 
(C.6b) 



both ingoing or both outgoing, and if both streams are 
relativistic relative to the no-going frame (as will be the 
case on the approach to the inner horizon) , then instead of 
the two expressions (|C.6p being approximately equal and 
of opposite sign, the two expressions are approximately 
equal and of the same sign. The ratio of the two equa- 
tions (jC.7p is then d\n f3/dlnua = — 1, and mass inflation 
is not ignited. 



Appendix D. Unexpected values of parameters 



Here a represents either of the ingoing or outgoing streams, 
and b represents the other stream. Inserting expressions (jC.5l) 
into equations (IB. 41) yields 
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-{-\ + P^ + Pbul) , (C.7a) 



2 (A - + y^bul) , (C.7b) 



where A is the same as before, equation pT|) . and the "ac- 
cretion rate" pb of each stream is given by equation 
generalizing equation (|T2|) . Note that it is important to 
distinguish the radii and rf, used as the independent 
variable on the left hand sides of equations (IC.7p . since 
the two streams follow different trajectories into the black 
hole. Equation (jC.7al) divided by equation (jC.7bp yields 
equation (|33p in the text. 



It is essential to the derivation of equations (IC.6P and 
(|C.7p from (|C.5p that one of the streams be ingoing and 
the other outgoing. If there are two streams that are 



It was commented at the end of §521 that the param- 
eters A, p,b, and pd of the analytic flt ([55]) to the model 
illustrated in Figure [3] are perhaps surprising given the 
boundary conditions of the model. The parameters are 
hstcd in Tabled 

One surprising feature is that the charge-to-mass ratio 
Q./Af, = 0.8 of the black hole set by the boundary con- 
ditions at the sonic point would predict A = 3 at the in- 
ner horizon of the corresponding Reissner-Nordstrom black 
hole, smaller than the value A = 14.6 actually measured 
in the model near the inner horizon. The reason for this 
is that the baryons are charged, and therefore repelled by 
the black hole, so the effective charge-to-mass of the black 
hole near the inner horizon is smaller than Q,/M, = 0.8, 
and A is consequently larger. 

A second surprising feature of the measured parameters 
is that the "accretion rate" parameter pd — 1-39 x 10^^ 
of the dark matter exceeds that pb — 1-93 x 10^'^ of the 
baryons, in spite of the fact that the boundary conditions 
at the sonic point set the dark matter density to be much 
less than the baryonic density, Pdj Pb — 10^'^. Ultimately, 
the reason for this is again that the baryons are charged 
not neutral. Although the dark-matter-to-baryon ratio 
Pdl Pb = 10^'^ at the sonic point may seem small, in fact 
it is almost as large as can be given the large charge-to- 
mass ratio Q,/M, = 0.8 of the black hole, and the small 
value w — 10^^ of the baryonic equation of state. If the 
ratio of neutral dark matter to baryons at the sonic point 
is increased, then the black hole is inclined to become 
more neutral, and it becomes impossible to achieve the 
desired charge-to-mass ratio Q./M, = 0.8 (as discussed by 
[49] , the mathematical condition that the parameters must 
satisfy is that the radial 4-gradient /3m must be spacelike 
at the sonic point outside the outer horizon). A small 
w = 10~^ makes it harder to achieve a large charge-to- 
mass ratio Q,/M, — 0.8, because then the sonic point 
is far from the black hole, and there is more room for 
the charged baryons to be repelled by the black hole as 
opposed to falling in. Thus one should think of the dark- 
matter-to-baryon ratio Pd/ Pb = 10"'^ of the model as being 
"large" rather than small. This is why the dark matter 
"accretion rate" pd measured near the inner horizon ex- 
ceeds the baryonic rate pb, Tabled] If the dark-matter-to- 
baryon ratio at the sonic point were set to a truly small 
value, then indeed the dark matter pd would be less than 
the baryonic pb near the inner horizon. 
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Appendix E. Mass accretion rate 

This paper parameterizes the accretion rate of the black 
hole as the rate M, of increase of its mass M, as mea- 
sured by distant observers, whereas parameterized 
the rate by a quantity 77^ = ar/t evaluated at the sonic 
point, where the boundary conditions are set. What is the 
relation between the two? 

As in [4^, the charge Q, and mass M, of the black 
hole at any instant are defined to be those that would be 
measured by a distant observer if there were no charge or 
mass outside the sonic point. 



), = Q and M, = M 



91 
2r 



at the sonic point . (E.l) 



The extra mass /2r added to the interior mass M is the 
mass-energy in the electric field outside the sonic point, 
given no charge outside the sonic point. 

In self-similar solutions, the black hole mass increases 
linearly with time, M, oc t, and the mass accretion rate 
M, is therefore 



M, = dM,/dt = M,/t , 



(E.2) 



where t is time measured at rest at infinity. 

The time t measured at infinity coincides with the 
proper time recorded on dust (dark matter) clocks that 
free-fall radially from rest at infinity. This can be seen 
as follows. First, for the line element ([5]), an interval 
of proper time recorded on a clock at rest in the tetrad 
frame is dt/a. This is true because a person at rest in 
the tetrad frame has, by definition, tetrad-frame 4- velocity 
u"* = {1,0,0,0}, so their coordinate-frame 4-velocity is 
uf' ^ em'^u'" = et^ = {a,/?t,0,0}, that is, dr/dt = Pt/a 
and dO/dt = d(f>/dt — 0, and it then follows from the line 
element ([5]) that the proper time interval of a person at rest 
in the tetrad frame is dt/a as asserted (the line element 
is in fact constructed to have this property). Second, 
for a freely-falling tetrad, the vierbein coefficient a is a 
function only of i, not r. This follows from equation (jA.7al) 
and the fact that, by definition, the proper acceleration 
ht experienced in a freely-falling frame vanishes. The 
fact that a{t) is a function only of t, not of r, expresses 
coordinate gauge freedom in the choice of time t. It is 
natural to fix the gauge by setting the coordinate time 
t equal to the proper time at rest at infinity. This is 
equivalent to setting a = 1 on dust clocks that free-fall 
from rest at infinity. In other words, as claimed, time t 
measured at rest at infinity coincides with the proper time 
recorded on dust clocks that free-fall radially from rest at 
infinity. To avoid any confusion, what this means is that if 
a dust clock falling through the sonic point says the time is 
t, and another dust clock falling through the sonic point a 
little later says the time is t+dt, then an observer at rest at 
infinity will say that the two instants are separated by time 
dt. One might worry about the difference in light travel 
times between sonic point and observer at the two instants. 



but for realistically small accretion rates the difference is 
negligible. 

In order to translate this time at infinity into what is 
happening at the sonic point, it is necessary to know the 
transformation between the dust frame and the baryonic 
frame at the sonic point, which is equivalent to knowing 
the vierbein coefficients /3„i of each. The coefficients for 
the baryonic frame are set by the boundary conditions at 
the sonic point. The interior mass at the sonic point, a 



gauge-invariant scalar, fixes f3f 



3'i in all frames. The 



only possible uncertainty is in the Pr of the dust at the 
sonic point. Technically, the dust /3r depends on the radial 
profile of mass-energy through which the dust has fallen. 
However, if the mass-energy outside the sonic point is 
neglected, an excellent approximation for realistically low 
accretion rates, then /3r = 1 at the sonic point, which is 
the value adopted here and in [i^, [l^ . 

By assumption, self-similar solutions possess conformal 
time translation invariance, which is to say that they re- 
main invariant under a scale transformation of t at fixed 
r/t. It follows that ^" defined by 



rrd. 



d 



dint 



d 



d_ 

dr 



{te"\ + re"V)9„ 



(E.3) 

is a dimensionless conformal Killing vector, or homothetic 
vector. Comparison to the inverse vierbein (|A.4p shows 
that the time component of the dimensionless homothetic 
vector is ^* — te^t/r — t/{ar). As argued above, the 
coefficient a equals one in the particular case of dust that 
free-falls from rest at infinity. Thus the time t at rest at 
infinity is 

* = <^ (E.4) 

where is the time component of the dimensionless ho- 
mothetic vector in the dust frame. Thus finally the mass 
accretion rate (IE. 21) is 



(E.5) 



M, = — r at the sonic point 



By comparison, the accretion rate parameter 77^ of (iol . [l6j 
1 



was 



Vs 



at the sonic point 



(E.6) 



where ^* is the time component of the dimensionless ho- 
mothetic vector in the baryonic frame. 
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